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This paper presents a heavy-traffic analysis of the behavior of a 
single-server queue under an Earliest-Deadline-First (EDF) schedul- 
ing policy in which customers have deadlines and are served only 
until their deadlines elapse. The performance of the system is mea- 
sured by the fraction of reneged work (the residual work lost due 
to elapsed deadlines) which is shown to be minimized by the EDF 
policy. The evolution of the lead time distribution of customers in 
queue is described by a measure-valued process. The heavy traffic 
limit of this (properly scaled) process is shown to be a deterministic 
function of the limit of the scaled workload process which, in turn, 
is identified to be a doubly reflected Brownian motion. This paper 
complements previous work by Doytchinov, Lehoczky and Shreve on 
the EDF discipline in which customers are served to completion even 
after their deadlines elapse. The fraction of reneged work in a heav- 
ily loaded system and the fraction of late work in the corresponding 
system without reneging are compared using explicit formulas based 
on the heavy traffic approximations. The formulas are validated by 
simulation results. 

1. Introduction. 
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1.1. Background and the reneging EDF model. In the last decade, at- 
tention has been paid to queueing systems in which customers have dead- 
hnes. Examples include telecommunication systems carrying digitized voice 
or video traffic, tracking systems and real-time control systems. In the case of 
voice or video, packetized information must be received, processed and dis- 
played within stringent timing bounds so that the integrity of the transmis- 
sion is maintained. Similarly, there are processing requirements for tracking 
systems that guarantee that a track can be successfully followed. Real-time 
control systems (e.g., those associated with modern avionics systems, man- 
ufacturing plants or automobiles) also gather data that must be processed 
within stringent timing requirements in order for the system to maintain sta- 
bility or react to changes in the operating environment. We refer to queueing 
systems that process tasks with deadlines as "real-time queueing systems." 

The performance of a real-time queueing system is measured by its abil- 
ity to meet the deadlines of the customers. This is in contrast to ordinary 
queueing systems in which the measure of performance is often customer 
delay, queue length or utilization of a service facility. We use the fraction 
of "reneged work," defined as the residual work not serviced due to elapsed 
deadlines, as our performance measure. To minimize this quantity, it is nec- 
essary to use a scheduling policy that takes deadlines into account. We use 
the Earliest-Deadline-First (EDF) policy, which reduces to the more familiar 
First-In-First-Out (FIFO) policy when all customers have the same dead- 
line. Under general assumptions, we prove that EDF is optimal with respect 
to this performance measure. A related result for G/M/c queues, in which 
the number of reneging customers is used as a performance measure, was 
obtained by Panwar and Towsley [29]. 

Heavy traffic analysis of a single real-time queue was initiated by Lehocz- 
ky [27]. This was put on a ffi'm mathematical foundation by Doytchinov, 
Lehoczky and Shreve (DLS) [7]. The accuracy of heavy traffic approxima- 
tions was developed in [22, 24]. The results of DLS were generalized to the 
case of acyclic networks in [23]. In these papers it was assumed that all 
customers are served to completion. The case in which late customers leave 
the system and their residual work is lost is addressed here. The main result 
of this paper is a heavy traffic convergence theorem, from which is derived 
a simple and practically useful approximation for the fraction of lost work 
when the system is heavily loaded. 

The mathematical formulation used by DLS and related papers is based 
on random measures. In addition to the usual queue length and workload 
processes associated with the queueing system, to model the evolution of a 
real-time queueing system, one must keep track of the lead time of each cus- 
tomer, that is, the time until the customer's deadline elapses. This is done by 
measure-valued queue length and workload processes. The measure-valued 
queue length process puts unit mass on the real line at the lead time of each 



EDF QUEUES WITH RENEGING 



3 



customer in the system, while the measure-valued workload process puts 
mass equal to the remaining service time of each customer at the lead time 
of that customer. These measures evolve dynamically as customers arrive, 
age and depart. Under the usual heavy traffic assumptions, since customers 
are served to completion in the DLS framework, it is easy to see that the 
ordinary scaled workload process converges weakly to a reflected Brown- 
ian motion with drift. DLS showed that the suitably scaled workload and 
queue length measure- valued processes converge to an explicit deterministic 
function of the workload process. 

In this paper customers leave the system when their deadlines elapse, 
which we refer to as reneging. Due to the preemptive nature of the EDF 
policy, it is not possible to determine at the time of admission whether a 
customer will be fully serviced before its deadline elapses. It is thus natu- 
ral to have the controller make the decision only at the time the deadline 
elapses. The system with reneging shows marked improvement in perfor- 
mance over the DLS system, in the sense that the fraction of reneged work 
in this system is much less than the fraction of work that becomes late in 
the DLS system. This improvement is because once a customer misses its 
deadline, the processor devotes no further effort to it, but rather turns its 
attention to customers that are not late. 

The system with reneging is considerably more difficult to analyze than 
the DLS system. In the reneging system, the evolution of the scalar total 
workload process depends on the entire lead time distribution of customers 
in queue and the nature of the EDF discipline. This is in stark contrast 
to the DLS system, where the total workload process is independent of the 
scheduling discipline, and is identical to that of any GI /G/1 queue with a 
work-conserving scheduling discipline. A key ingredient of our analysis is a 
mapping on the space of measure-valued functions which, when applied to 
the DLS system, yields another system (that we call the reference system) 
whose difference from the reneging system vanishes in heavy traffic. This 
mapping can be viewed as a generalization of the scalar double reflection 
map to measure-valued processes, and, using its continuity properties, we 
identify the heavy traffic limit of the reference and hence the reneging sys- 
tems. Specifically, we show that the limit of the scaled workload process is a 
doubly reflected Brownian motion with lower barrier zero and upper barrier 
at the mean of the lead time distribution. We also show that, conditional 
on the limiting workload, the resulting limiting measure-valued workload 
process is the same limiting process as when customers are served to com- 
pletion, that is, in the DLS system. However, the workload processes in 
these two systems differ, and so the unconditional limiting lead-time profiles 
of these two systems differ accordingly. In particular, unlike in the DLS sys- 
tem, the measure-valued workload process in the reneging system is always 
concentrated on the positive real line due to the absence of late work in the 
reneging system. 
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1.2. Prediction formulas. The results of this paper suggest a simple for- 
mula for the fraction of lost work in the EDF system with reneging. In 
particular, consider a single-server queue with traffic intensity p = A//i that 
is near one, where 1/A is the mean interarrival time and l//i is the mean 
service time. Let a and /3 be the standard deviations of the interarrival times 
and service times, respectively, and set o"^ = A(a^ + /3^), which we assume is 
nonzero. Let D denote the mean lead time for arriving customers. Finally, 
set 9 = 2(1 — p)lo^ . Under these circumstances, 



(1.1) Fraction of lost work in reneging system ~ e 



l-p 



p{l - e 



This formula is derived in Section 7.1 and compared with simulations in 
Section 7.2. If /3= 1, in place of (1.1) we have 

(1.2) Fraction of lost work in reneging system ~ 

Analysis of the limit of the standard (nonreneging) system suggests that 
when p<l [see (7.8) and (7.9)], 

(1.3) Fraction of late work in standard system w e~^^, 
which, together with (1.1), yields the approximation 

^ ^ Lost work in reneging system 1 ~ P 

Late work in standard system p{\ _ e~^^) 

If /9 > 1 then all work is late in the limiting standard (nonreneging) system, 
which leads to the approximation 

Lost work in reneging system 
Late work in standard system 2D 
When plotted on a log scale, the fraction of lost work in the reneging system 
and the fraction of late work in the standard system will be linear in D, 
provided that e^^ ^ 1, and these two plots will be separated by log((l — 
p)/p)- When performance is measured in terms of the work whose service 
requirement is not met by the time its deadline elapses, then the reneging 
system is far superior to the nonreneging system. We refer the reader to the 
simulations in Section 7.2. 

The situation with reneged customers as opposed to reneged work is more 
complicated. DLS shows that the number of customers in the limiting stan- 
dard system at any time is just A times the amount of work, the number of 
late customers is A times the amount of late work and hence 

Fraction of late customers in the standard system 

(1.6) 

~ Fraction of late work in the standard system 
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[see also (7.8) and its derivation for the case p <1]. In the Umiting reneg- 
ing system, the number of customers who arrive by a certain time and the 
number of customers in system at that time is A times the amount of ar- 
rived work and A times the amount of work in the system (Corollary 3.7), 
respectively, but the number of customers who renege by a certain time is 
not necessarily A times the amount of reneged work by that time (see Re- 
mark 7.2). In particular, we do not have a formula like (1.6) for the reneging 
system. If the arrival process is Poisson, the fraction of lost customers in the 
reneging system can be estimated by a heuristic argument [see (7.7)] which 
gives instead 

Fraction of lost customers in reneging system 

(1.7) 

2 

^ ~o7y) ^ (Fraction of lost work in reneging system). 

/x^p^ -|- 1 

1.3. Related work and outline of paper. Measure-valued processes have 
recently gained prominence in queueing theory. Decreusefond and Moyal 
[5] use such processes to obtain the fluid limit of an EDF M/M/1 queue 
with reneging. Unlike our scaling (2.4) of lead times by ^/n, they scale lead 
times by n and obtain a characterization of the limiting lead-time measure- 
valued process via a transport equation. In a different setting. Ward and 
Glynn [33, 34] find limits of FIFO queues with reneging. Measure- valued 
processes have also proved useful in the heavy traffic analysis of queues 
with scheduling disciplines other than EDF such as last-in- first-out [28], 
processor sharing [11, 12], and shortest remaining processing time [6, 13]. 
As dynamical systems, queueing systems present a mathematical challenge 
due to discontinuities in their evolution at boundaries (which denote empty 
queues). The heavy traffic analysis of queueing systems described by M""- 
valued processes has been facilitated by the use of representations in terms 
of continuous mappings on M" [4, 8, 14, 31, 36]. This work demonstrates that 
this perspective can also be useful when the queueing system is represented 
by a more complicated, measure-valued process (see also [18] for recent work 
that takes a similar perspective). 

Section 2 introduces our model. Section 3 summarizes the main results, 
and proofs of these results are given in Section 6. Section 4 introduces the ref- 
erence workload process and its decomposition, and describes its evolution. 
This reference workload process is easier to analyze than the workload pro- 
cess with reneging but the two are shown to have the same asymptotic behav- 
ior. Comparisons between the reference workload process and the reneging 
workload process are presented in Section 5. Section 7 presents simulation 
results. A proof of optimality of EDF, that may be of independent interest, 
is in the Appendix. 
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2. The model, assumptions and notation. 

2.1. Notation. Let M be the set of real numbers. For a, 6 G M, a V 6 is the 
maximum of a and b, a A 6 is the minimum and is the maximum of a and 
0. Also, inf{0} should be understood as +oo, while sup{0} and max{0} 
should be understood as — oo. Moreover, if a < 6, then the interval [b,a] is 
understood to be 0. 

Denote by M the set of all finite, nonnegative measures on ,S(M), the Borel 
subsets of M. Under the weak topology, is a Polish space. We denote the 
measure in A4 that puts one unit of mass at the point rr G M, that is, the 
Dirac measure at x, by 5x- When v £ Ai and B is an interval {a,b] or a 
singleton {a}, we will simply write i/(a, 6] and z/{a} instead of i^((a, &]) and 

Let T > be given. Given a Polish space X, we use Dx[0,(x>) (resp., 
Dx[0,T]) to denote the space of right-continuous functions with left-hand 
limits (ROLL functions) from [0,oo) (resp., [0,T]) to X, equipped with the 
Skorokhod Ji topology. See [9] for details. When dealing with Dx[0,oo) or 
Dx[0,T], we typically consider X = M or M^, with appropriate dimension 
d for vector- valued functions, or X = Ai, unless explicitly stated otherwise. 
When X = R or M, for t > and x G Dx[0, oo), we write x{t-) for the left- 
hand limit limsi-t x{s), and we define Ax{t) to be the jump in x at time t, that 

is, Ax{t) = x{t) — x(t—). Finally, given Z)x[0,oo)- valued random variables 
Zn,n £ N, defined, respectively, on the probability spaces {Qn,J^n,^n), n G 
N, and a Dx [0, oo)-valued random variable Z defined on a probability space 
J^, P), we say Z^"^^ converges in distribution to Z and write Z^ =^ Z, if 
for every bounded continuous function / on Dx[0,oo), lim„_5.oo E„[/(Z„)] = 
E[/(Z)]. Here E„ and E are expectations taken with respect to and P, 
respectively. 

2.2. The model with reneging. We have a sequence of single-station queue- 
ing systems, each serving one class of customers. The queueing systems 
are indexed by superscript (n). The inter-arrival times for the customers 

are a sequence of strictly positive, independent, identically dis- 

tributed random variables with common mean t4-y and standard deviation 
Q,{"). The service times are {^^"^Ij^i, another sequence of positive, indepen- 
dent, identically distributed random variables with common mean and 

standard deviation 

If the initial condition of the nth queue were not zero, then we would need 
to specify an initial workload measure-valued process and frontier [these 
terms are defined in (2.17) and (2.19) below] in such a way that these have 
limits under the heavy traffic scaling. However, if the limit of the initial 
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scaled workload process were not of the form appearing in Theorem 3.2 
below, then the workload process would be expected to have a jump at time 
zero. To avoid these complications, we assume that each queue is empty at 
time zero. 

We define the customer arrival times 

(2.1) S(")^0, 5(")^^n;"\ k>l, 

i=l 

the customer arrival process 

(2.2) ^(")(t)=max{A;;5^"^<t}, t>0, 
and the work arrival process 

[t] 

(2.3) y{«)(t)^^^;W, t>0. 

The work that has arrived to the queue by time t is then y(")(A(")(t)). 

(n) 

Each customer arrives with an initial lead time Lj , the time between 
the arrival time and the deadline for completion of service for that customer. 
These initial lead times are independent and identically distributed with 

(2.4) ¥{L^;^ <V^y} = G{y), 

where G is a right-continuous cumulative distribution function. We define 

(2.5) y*=inf{yGM|G(y)>0}, = min{y G M|G(y) = 1} 

and assume that < y* < y* < +oo. We assume that for every n, the se- 
quences {Uj^^}^i, {v^"^}j^i and {L^J^^}^^^ are mutually independent. See 
Remark 3.9 for a discussion of these assumptions. 

We assume that customers are served using the Earliest-Deadline-First 
(EDF) queue discipline, that is, the customer with the shortest lead time 
receives service. Preemption occurs when a customer more urgent than the 
customer in service arrives (we assume preempt-resume). There is no set up, 
switch-over, or other type of overhead. If the j'th customer is still present 
in the system (either waiting for service or receiving it) when his deadline 

(n) (n) 

passes, that is, at the time Sj + Lj , he leaves the queue immediately. 
This may be interpreted as either reneging or the result of an action of an 
external controller. 

We define W^"\t), the workload process at time t, as the remaining pro- 
cessing time of all the customers in the system at this time. We define Ry^ {t) 
to be the amount of work that reneges in the time interval [0, t] . The queue 
length process Q^^\t) is the number of customers in the queue at time t. 
The queueing system described above will be referred to as the EDF system 
with reneging. 
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2.3. The standard EDF model. We also have a sequence, indexed by su- 
perscript (n) , of standard EDF systems, with the same stochastic primitives 
as the EDF systems with reneging. In each of these standard systems, the 
server serves the customer with the shortest lead time, preemption occurs 
as in the reneging system, but late customers (customers with negative lead 
times) stay in the system until served to completion. The performance pro- 
cesses associated with the standard system will be denoted by the same 
symbols as their counterparts from the system with reneging, but with ad- 
ditional subscript S. For example, Wg^\t) denotes the workload in the stan- 
dard system at time t. The arrival processes and are the same 
for the both systems, so we will not attach the subscript S to them. 

The standard EDF system is easier to analyze than the EDF system with 

(n) 

reneging in several ways. For instance, the workload Wg in the standard 
system coincides with the workload of a corresponding G/G/1 queue (with 
the same primitives) under any nonidling scheduling policy. More precisely, 
in the standard system the netput process 

(2.6) iv(")(t) = y(")(^(")(t))-t 

measures the amount of work in queue at time t provided that the server is 
never idle up to time t, and the cumulative idleness process 

(2.7) 4")(i)^- inf iV(")(s) 

0<s<t 

gives the amount of time the server is idle. Adding these two processes 
together, we obtain the workload process for the standard system 

(2.8) VT^"^ (t) = (t) + /J') (t) . 

(All the above processes are RCLL.) In contrast, the evolution of the work- 
load Ty(") in the reneging system is more complex and depends not only 
on the residual service times but also on the lead times of all customers in 
the queue. Our analysis of the reneging system will be facilitated by results 
from [7] on the heavy traffic analysis of the standard EDF system. 

2.4. Heavy traffic assumptions. We assume that the following limits ex- 
ist: 

lim A(") = A, lim = A, 

71— ^-OO 71— >00 

(2.9) 

lim = Q, lim = f3, 

n— >oo 71— ^-oo 

and, moreover, A > and + /3^ > 0. Define the traffic intensity p^") = 
We make the heavy traffic assumption 

(2.10) lim V^(l-/0^"^)=7 
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for some 7 G M. We also impose the Lindeberg condition on the inter-arrival 
and service times: for every c > 0, 



limE[(.4">-(A("))-')^I,,„,.,_,,,,,_.,^^^, 



(2.11) 



71— >00 



„iT„E[(„;»'-(,<»))-'fi„,„_,^„,„_.,^^^,i=o. 



We introduce the heavy traffic scaling for the idleness process in the stan- 
dard system and the workload and queue length processes for both EDF 
systems 



\/n -v/n 



Q('^)(t) = ^g('^)(nt) 
\/n 



and the centered heavy traffic scaling for the arrival processes 



The scaled netput process (which is the same for both systems) is given by 
(2.12) (t ) = ^ ( A^'^) (nt) )-nt]. 



n 

Note that, by (2.8), W^''\t) = N^^-^t) +I^p\t). 

It follows from Theorem 3.1 in [30] and Theorem 7.3.2 in [36] that 

(2.13) (5("),1("))^(5*,A*), 

where A* is a zero-drift Brownian motion with variance a^X^ per unit time 
and 

(2.14) S*iXt) = -\A*it), t>0. 

A 

It is a standard result [16] that 

(2.15) (7V("),4"\t?("))^(iV*,/^,W^^), 
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where A^* is a Brownian motion with variance (q^ + /3^)A per unit time and 
drift —7, 

(2.16) Psit) ^- mm N*{s), W*s{t) = N* (t) + 1^- 

0<s<t 

In other words, Wg is a Brownian motion reflected at with variance (a^ + 
/3^)A per unit time and drift —7 and Ig causes the reflection. 

2.5. Measure-valued processes and frontiers. To study whether tasks or 
customers meet their timing requirements, one must keep track of customer 
lead times. The action of the EDF disciphne requires knowledge of the cur- 
rent lead times of all customers in system. We represent this information via 
a cohection of measure-valued stochastic processes. 

Customer arrival measure-valued process: 

A^^'^ (t)(B) — / arrivals by time t, whether or not still in the 

\ system at time t, having lead times at time t in S € B(M 

Workload arrival measure-valued process: 

V^^Ut)(B) — J^*-"^^ arrived by time t, whether or not still in the sys- 
I tem at time t, having lead times at time t in S S B{M.) 

Queue length measure-valued process: 

Q^^^t)(B) — I ^^^^^^ customers in the queue at time 1 
^ \^t having lead times at time t in B G B{W) J ' 

Workload measure-valued process: 



(2.17)W(")(t)(5) = <^ 



A J Work in the queue at time t associated with cus- 
tomers having lead times at time t in B £ B(M) 



The latter two processes describe the behavior of the EDF system with 
reneging. Their counterparts for the standard EDF system will be denoted by 
Q^g\t) and >V^"'^(i), respectively. The following relationships easily follow: 

A('^)(t) = ^(")(t)(M), F(")(A(")(t)) = v(")(t)(]R), 

Vr(")(f) = W(")(t)(0,oo), Q(")(t) = Q(")(t)(0,oo), 

M^(")(t) = W^")(t)(M), Q^^\t) = Q^^\tm. 

In addition, we can represent the reneged work as follows: 

(2.18) R^^\t)= w(")(.-){0}. 

0<s<t 
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In order to study the behavior of the EDF queue discipUne, it is useful 
to keep track of the largest lead time of all customers, whether present or 
departed, who have ever been in service. We define the frontier 

{The maximum of the largest lead time of 
all customers who have ever been in service, 
whether still present or not, and y/ny* — t 

for the EDF system with reneging, and its counterpart F^"^ (t) for the stan- 
dard EDF system. Prior to arrival of the first customer, F^"'\t) and F^'^\t) 
equal y/ny* — t. For the EDF system with reneging, we define the current 
lead time 

— f I^^^d time of the customer in service 
I or F^") (t) if the queue is empty 

In the reneging system, there is no customer with lead time smaller than 
C^"^(t), and there has never been a customer in service whose lead time, if 
the customer were still present, would exceed F^'^\t). Furthermore, C^"^(t) < 

for all t > 0. The processes C^''\ F('^) and F^ are RCLL. 
We introduce heavy traffic scalings. For the real-valued processes 
CW,F("),F_5"\^('^\Q('^),i2[;^\ we define = ^Z(")(nt) and for the 

measure- valued processes = Q^*") , W^") , qJT^ , , , V^") , we define 
Z'-'^\t){B) = -^Z('^\nt){^B) for every Borel set B CR. 

3. Main results. Before stating our main results, we summarize the re- 
sults for the standard EDF system that were obtained in [7] — in particular, 
we recall Proposition 3.10 and Theorem 3.1 of [7] which characterize the lim- 
iting distributions of the workload measure and the queue length measure 
in the standard system. Let 

(3.1) H{y)^ Til- Gir,))drj^^ I ^ ^ 




if y >y*. 

The function H maps (— oo,y*] onto [0,oo) and is strictly decreasing and 
Lipschitz continuous with Lipschitz constant 1 on (— oo,?/*]. Therefore, there 
exists a continuous inverse function that maps [0,oo) onto (— oo,?/*]. 

Proposition 3 . 1 (Proposition 3. 10 [7] ) . We have F^"^ =^ F| as n ^ oo, 
where the limiting scaled frontier process Fg for the standard EDF system 
is explicitly given by 

(3.2) F*sit) = H-\wm, t>0, 
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with Wg equal to Brownian motion with variance {a^ + /3^)A per unit time 
and drift —7, reflected at 0. 

Theorem 3.2 (Theorem 3.1 [7]). Let and Q*g be the measure-valued 
processes defined, respectively, by 

(3.3) W*sit)iB)= [ il-Giy))dy, Q^tm ^ XW*sit)iB), 

JBn[F*{t),oo) 

for all Borel sets 5 C R. Then wf' and ^Q%, as 00. 

Remark 3.3. The proofs in [7] can be modified to show that the con- 
vergences in (3.3) are in fact joint, that is, {W^^\Q.P) (W^, Q*s). 

There is lateness in the standard EDF system if and only if the measure- 
valued workload process has positive mass on the negative half line. Theo- 
rem 3.2 shows that, in the heavy traffic limit, this occurs exactly when the 
limiting scaled frontier process lies to the left of or, equivalently (by 

Proposition 3.1), when is greater than H{0) = E[L^"'''/-y/n], the mean of 
the scaled lead-time distribution. In the reneging system, there is no late- 
ness, and the amount of work that reneges is precisely the amount required 
to prevent lateness. Thus it is natural to expect that the limiting workload 
in the reneging system will be constrained to remain below H(0). Let W* 
be a Brownian motion with variance (a^ -|- /3^)A per unit time and drift —7, 
reflected at and H{0). The first main result of this paper is that W* is the 
limiting workload in the reneging system. 

Theorem 3.4. As n — >■ 00, 

The next two results of this paper are the following counterparts of Propo- 
sition 3.1 and Theorem 3.2 for the EDF system with reneging. 

Proposition 3.5. We have F^^^ ^ F* as 00, where 

(3.4) F*{t) = H-\W*{t)), t>0. 

In other words, the process F* defined by (3.4) is the limiting scaled 
frontier process for the EDF system with reneging. 

Theorem 3.6. Let W* and Q* be the measure-valued processes defined 

by 

(3.5) W*{t){B)^ f {l-G{y))dy, Q* {t){B) ^ \W* {t){B), 

JBn[F*{t),oo) 

for all Borel sets 5 CM. Then (W("),Q(")) (W*,Q*) as n-^ 00. 



EDF QUEUES WITH RENEGING 



13 



By Theorem 3.6, the total masses of >V(") and Q(") must converge jointly 
to the total masses of W* and Q'^'^\ respectively. Substituting i? = M in (3.5) 
and using (3.1) and (3.4), we see that W*(i)(M) =H{F*{t)) = W*{t) and 
we recover Theorem 3.4. In fact, we have a stronger result. 

Corollary 3.7. ^sn^oo, (!?("), Q(")) ^ (Vl^*,AT^*). 

Theorem 3.6 also shows that the limiting instantaneous lead-time profiles 
of customers in the EDF system with reneging conditioned on the value of the 
(limiting) workload in the system are the same as in the case of the standard 
EDF system. However, the limiting real-valued workload process for the 
EDF system with reneging is W* , the doubly reflected Brownian motion 
and the unconditional limiting lead-time profiles for these two systems differ 
accordingly. 

We also have a characterization of the limiting amount of reneged work. 

Theorem 3.8. As n^oo, =^R^, where is the local time at 
H{0) of the doubly reflected Brownian motion W* . 

Although these results are intuitive in light of the behavior of the standard 
EDF system, the proofs are challenging. Moreover, counter to what one 
might expect, the result for queue lengths analogous to Theorem 3.8 is false. 
Specifically, although Corollary 3.7 shows that Q^^^ converges to the doubly 

reflected Brownian motion Q* = XW* on [0, XH{0)], the scaled sequence 

^(n) 

Rq , n G N, of reneged customers does not converge to the local time Ai?^ 
of Q* at XH{0). This observation, which is elaborated upon in Section 7, 
emphasizes the need for a rigorous justification of intuitive statements. 

The proof of Theorem 3.4 is in Section 6.1.1, the proofs of Proposition 3.5 
and Theorem 3.6 are in Section 6.1.2, and Section 6.2 contains the proof of 
Theorem 3.8. We also establish an optimality property for EDF, Theorem 
5.1. 

Remark 3.9. The assumption made in (2.5) that the support of the 
lead time distribution is bounded above by y* < oo is mainly technical. It 
is expected that the analysis in [21] for the standard EDF system under a 
weaker second moment condition can be applied to the reneging system as 
well. On the other hand, the lower bound y^, > on the lead time distri- 
bution or some restriction on the behavior of the density of the lead time 
distribution at appears to be necessary. Indeed, the work of Ward and 
Glynn [33, 34] on FIFO queues with reneging suggests that in the absence 
of such an assumption, the limiting workload process may no longer be a 
reflected Brownian motion, and its properties may exhibit strong sensitivity 
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to the density of the lead-time distribution near 0. From a modehng point 
of view, it is reasonable to impose a strictly positive lower bound > so 
as to avoid nonnegligible "intrinsic lateness," in which an arriving customer 
has such a small initial lead time that he would be late even if there were 
no other customers in the system. 

In [21] the assumption of independence between the sequence of interar- 
rival times and lead times is also removed and a more complex version of 
Theorem 3.2 is obtained. Starting from that more complex result, the limit 
of the reneging system can be obtained along the lines of this paper. 

4. The reference system. In this section we introduce an auxiliary ref- 
erence workload measure-valued process and the corresponding real- 
valued reference workload process U^"^^ . In the special case of constant ini- 
tial lead times (i.e., =y*), in which EDF reduces to the well-known FIFO 
service discipline, ZY^"^ and ?7^"^ coincide with W^"^ and W^'^\ respectively. 
In general, these processes do not coincide (see Example 4.6), but, as we will 
show in Section 6.1, the difference between the diffusion-scaled versions of 
and W^""^ is negligible under heavy-traffic conditions. The advantage of 
working with the reference system, rather than the reneging system, is that 
can be represented explicitly as a certain mapping $ of the measure- 
valued workload process W^"^ in the standard system. As shown in Section 
6.1, continuity properties of the mapping $ enable an easy characterization 
of the limiting distributions of ZY*-"^ and in heavy traffic. 

We begin with Section 4.1, where we define the reference system and 
provide a useful decomposition of the process C/W. In Section 4.2 we provide 
a detailed description of the evolution ofU^^h 

4.1. Definition and properties of the reference workload. In Section 4.1.1, 
we introduce a deterministic mapping on the space of measure- valued func- 
tions that is used to define the reference workload. Then, in Section 4.1.2, 
we provide a decomposition of the reference workload process. 

4.1.1. A mapping $ of measure-valued processes. We define a sequence 
of reference workload measure-valued processes for the EDF system with 
reneging by the formula 



(4.1) =$(WJ'^), 
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for every G D_a4[0,oo), t>0 and y G M. (The claim that ^ does indeed 
map D_m[0,oo) into D_\4[0,oo) is justified in Lemma 4.1 below.) We also 
define the (real-valued) reference workload process C/^"^ as the total mass of 
that is, 

(4.3) C/(")(t)=ZY(")(t)(]R) VtG[0,oo). 

(n) 

The frontier Fg defined in Section 2.3 played a crucial role in the descrip- 
tion and analysis of the evolution of the standard system in [7] . In a similar 
fashion, it will be useful to define the reference frontier 

(^A^ p{n).,^ A rinf{yG]R|ZYW(t)(-oo,y]>0}, if C/(")(i) > 0, 
i4.4j ^*^-\+oo, ifC/H(t) = 0. 

By definition, (t) is the leftmost point of support of the random measure 
U^'^\t) [understood as oo if Z^(")(t) = 0]. The process -E^") has ROLL paths. 
From (4.1)-(4.3) we have 

(4.5) U^-\t){-oc,y] = [w(")(t)(-oo,y] -i^(")(t)]+, 

(4.6) (t) = Vf^^"^ {t) - (") {t) , 
where 

(4.7) i^(")(t)= max(w5'^(s)(-oo,0] A inf wi''\u)}. 

se[o,i] I ue[s,t] ^ ) 

In (4.7) we may write maximum rather than supremum because the pro- 
cess w'g'\-){—oo,0] never jumps down. Note from (4.6) and (4.7) that 
< K(")(t) < W^''\t) and so for all t > 0, 

(4.8) (t)<VF5(t) • 

According to (4.6), the reference workload process C/^") is the standard work- 
load process W^"^ with mass K^"^^ removed. Equation (4.5) shows that this 

(n) 

mass is removed from the left-hand side of the support of Wg . Moreover, 
since U^"'\t){—oo,y] > for all y to the right of the frontier E^"'\t), it is 
clear from (4.1) and (4.2) that for t £ [0,oo), 2/2 > yi > E^''^\t), 

(4.9) ZY(")(t)(yi,y2]=Z^(")(t)(-oo,y2]-^/("^(t)(-oo,yi]=w(")(t)(yi,y2], 

which shows that Z//^"^ coincides with W^"^ strictly to the right of E^"'\ 

In the following lemma, we establish some basic properties of ^> that show, 
in particular, that U^"'\t), t>0, and U^"'\t), t>0, are stochastic processes 
with sample paths in D_a4[0,oo) and Z)ir^[0, 00), respectively. Although ^> 
is not continuous on D_a4[0,oo), the lemma shows that it satisfies a certain 
continuity property that will be sufficient for our purposes. 
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Lemma 4.1. For every t G [0,oo), <I>(^)(i)(— oo, 0] = 0. Moreover, <I> maps 
Dm [0, oo) to Dm [0, oo) . Furthermore, if a sequence n G N, in Dj^ [0, oo) 
converges to fj,G D_\4[0,oo), where fi is continuous and for every t G [0,oo), 
fj,{t){0} =0, then ^{Hn) converges to in D_m[0,oo). 



Proof. The first statement follows from the simple observation that, 
due to the nonnegativity of /.i and (4.2), 



< $(/i)(t)(-oo, 0] < [fi{t){-oo, 0] - n{t){-oo, 0] A /x(t)(M)]" 



0. 



Also, since the right-hand side of (4.2) is nondecreasing and right-continuous 
in y, we know that ^>(/i)(t) G 7W for every t > 0. Now, observe that <I>(/i)(t) = 

^{fi{t),T{fi){t)), where ^ : x R is the mapping ^'(i/, j;)(-oo, y] = 

(i^(— oo, y] — x)"*" for all y G M and F : Dj^ [0, oo) i— M is defined by 



r(/_i)(t) = sup (/_i(s)(— 00,0] A inf fi{u 

s6[0,tp ue[s,t] 



Vt G [0,oo). 



Using the fact that weak convergence of measures on M is equivalent to 
convergence of the cumulative distribution functions at continuity points 
of the limit, one can verify that ^ is continuous on 7W x M. To show that 
$(/i) G D;k[0,oo), it suffices to show that T{fi) G D[0,oo). For this, we fix 
t G [0, oo) and write 



r(M)(t + e)-r(^)(i) 



sup 



/i(s)(— OO, 0] A inf fj.{u 

ue[s,t] 



A inf fi{u){R) \J Z{ii,e){t) 

w6[i,t+£] 



— sup 

se[o,t] 

where we define 



/i(s)(— oo, 0] A inf /i(u 

«e[s,t] 



Z{ii,e){t) 



sup 

se[i,t+( 



/i(s)(-oo,0] A inf ^(u)(M) 



Since /i G D;\4[0,oo) implies ^{u) converges weakly to /i(t) as u|t, we have 
lim„j_t //(«)(R) = ;u(i)(]R) and /i(t)(— oo,0] > limsup^^j ^(s)(— oo, 0] by Port- 
manteau's theorem. This, in turn, implies that lim£_j.o Z{^, e)(i) = Ai(i)(— co, 0] 
for all t > 0. Combining the above properties, it is easy to deduce that 
r(/x)(t + e) — r(//)(t) — >• as e I 0, and the right-continuity of $(/i) follows. 
The existence of left limits for r(/u), and hence for ^{u), can be established 
by an analogous but simpler argument. 

Now, suppose ^ri converges to /i in D_a4[0,oo) and /u is continuous with 
^(t){0} = for every t > 0. Then finit) converges weakly to /i(t) uniformly 
for t in compact sets (u.o.c.) (see [2]). Since is a continuity point for /x(t). 
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this implies oo,0] and /i„(t)(M) converge u.o.c. to oo,0] and 

respectively. This shows that r{ij,n){t) converges u.o.c. to T{ii){t), 
which, when combined with the continuity of ^, shows that <^(;U„)(t) con- 
verges weakly u.o.c. to In particular, this shows converges to 
<^{l^i) hi Dm[0,oo). □ 

As an immediate consequence of the lemma, the definitions of ZY^"^ and 
and the fact that is a purely atomic measure, we have, for all 

t>0, 

(4.10) ZY(")(i)(-oo,0] = and E^''\t)>0. 

4.1.2. A decomposition of the reference workload. We establish a decom- 
position of K^"^^ into its increasing and decreasing parts. Define cJq"^ = and 
W^^'^ (0—) = 0. For A; = 0, 1, 2, ... , define recursively 

) ^ min{t > v max W^^\s){-oo,0] 

(4.11) 



(4.12) ^ min{t > rt^\W^-\t) > W^-\t-)}. 

In addition, for t G [0,oo), define 

K^-\t) = (4"^-) V max w(-)(.)(-oo,0] 

(4.13) 

-h4'^^(4"^-)_, 
K^-\t) = - E[(^"^(^S) - (4"^ ^ * - -i-i))^ 

fceN 

(4.14) 

-W4")(ra)]. 

Theorem 4.2. We have 
(4.15) K^'^) = ^ - 

where K^'^ and K^'^ are the positive and negative variations of K^"'\ More- 
over, 



(4.16) / I^^(„)(,)>o}'^^^Hs) = 0. 

[0,oo) 
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The theorem is easily deduced from Propositions 4.3 and 4.4 and Remark 
4.5 below. The rest of the section is devoted to establishing these results. 

Observe that the late work Wg^\s)(—oo,0] is right-continuous in s, re- 
maining constant or moving down at rate one and jumping up. There- 
fore, the maximum on the right-hand side of (4.11) is obtained. Addition- 

(n) (n) 

ally, because of the right-continuity of Wg and Wg , the minimum in 

this equation is also obtained. Finally, W^"''*(s)(— oo, 0] can never exceed 

Wg^\s) = yv'g'\s)(M.), and W^^^ never jumps down, so we must in fact 
have 

(4.17) m4")(4")-)V max w(")(.)(-oo,0]=Ty^")(rf)). 

For k>l, is the first arrival time after t^"\ . We thus have 

(4.18) wt\t) = {W^^\rtJ,)-{t-rtJ,))\ rtJ.KtKai^K 
We further have 

(4.19) = ai") = ri") < aj") < ^ < < • • • . 



Proposition 4.3. For each k>l, we have 

(4.20) if(")(t) = T^^"^(4"^-) V max wL"^(s)(-cx),0] 

for t G [<7^"^ , T^^] . In particular, K^"^^ is nondecreasing on the interval [<t^"^ , 7"^"^] • 
Proof. We proceed by induction on k. For the base case k = 1, note 

(n) 

that the standard EDF system is empty before the time a\ . Therefore, 
Wg^\a^^ — ) = 0, and to prove (4.20), we must show that 

(4.21) i^(")(t)= max W^"^(s)( -00,0], cj^") < t < rf \ 

For t G [cjj"\ r}"^], let s^''^\t) be the largest number in [cr^"\t] satisfying 

(4.22) Wj'^(s(")(t))(-oo,0] = max W^''^(s)(-oo, 0]. 

se[o,t] 

For u G [s^"\t),t], we have 

wP{s^''\t)){-oo,0]= max W^"^(s)(-oo,0], 
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which is less than or equal to VFg"^ (u) by the definition of t[^^ and equation 
(4.17). Therefore, 

maxM"^(s)(-oo,0] =wL"^(s(")(t))(-oo,0] < inf W^''\u). 

s&[0,t] wG[s(")(t),t] 

Equation (4.21) follows from (4.7). 

We assume (4.20) holds for some k and prove it for k + 1. For t £ [(t^^^i, Tj,"^-^] 

K^"-\t)= max |w^"^(s)(-oo,0] A inf wi"'\u)} 
^£[0,4';',)'^ ^ J 

(4.23) 

V max |wi"^(s)(-oo,0] A inf wi''\u)]. 
Equation (4.20) with k replaced by A; + 1 will follow once we show that 



(4.24) max f wl"Vs)(-oo, 0] A inf wi''\u)} 



max |w^"Vs)(-oo,Ol A inf VF^"Vu)| 



and 
(4.25) 

= max W^"^(s)(-oo,0]. 

For (4.24), we observe that because W^"^(s)(— oo,0] and infs<„<f Wg^\u) 
regarded as functions of s, cannot increase except by a jump, the maximum 
the : 
An) 



on the left-hand side of (4.24) is attained. Let s^"^ be the largest number in 



[0, o"^._i_-i^) attaining this maximum. We have 
max |w^"Vs)(-oo,0] A inf wi''\u)\ 



se[o 



:W(-)(4"))(-oo,0]A irrf W^-\u) < W^-\u) e [s^-\ a^i,) 

ne[4"',i] 



and so 

(4.26) max M^\s){-oo,0] A inf Tvi"^(n)| < VF^"^(a["^ 



On the other hand, by the inequalities r^"^ < o"^'^^-^ <t< t^"\, definition 
(4.7), the induction hypothesis, and equation (4.17), we have 

max < 



max |wl"Vs)(-oo,0l A inf wlp\u)} 
-JnJ^U^ ue[s,t] ) 
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> max |wi"^(s)(-oo,0] A inf W^!'\u) A inf W^''\u)} 

= i^W(rf))A inf W^-\u) 

(wi''\ai'^-)V max W^"^(s)(-oo, 0]) A inf T^J"^(n) 



l^^")(rf))A inf VF|"Hn)= inf 



Equation (4.18) implies Wf>{u) > W^""' {a]^_l^-) for r^"^ <u< For 
c^i+i < ti < t < r^"\, (4.11) implies that 

W^5"^(4+\-)V max wP{s){-^,0]<W^''\u), 

and so again we have Wg"'\u) > ^^^"■'((T^'^^-^— ). Finally, if u = t = t^"\, then 
(4.17) implies that > It follows that 

inf W^-\u)>W^-\a^-l,-). 

This gives the reverse of the inequality (4.26), and thus (4.24) is proved. 
For (4.25), we let t^"^ attain the maximum in max^^j^(„) >V_5"^(s)(— oo,0]. 

For u G [t^^^t], we have from (4.11) and (4.17) that 

Wi")(4"))(-oo,0]= max {s){-oc,0] <W^''\u), 

se[4+i,-"] 

and hence W^"^(4"^)(-oo, 0] < mf^g[j{n) It follows that 

max |wL"^(s)(-oo,0] A inf wl"^(u)| 

<n'f (4"^)(-oo,0]=w(")(4"))(-oo,0]A inf T^'^)(n) 

< max |wj'^(s)(-oo,0] A inf wi''\u)}, 

which establishes (4.25). □ 

Proposition 4.4. For each k > 1, we have 
(4.27) K(-\t) = {W^-\ri-J,)-it-rj::^,))\ ri"i<t<a^^\ 
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In particular, K^") is nonincreasing on [t^."\ , c^."'' ) . 

Proof. For all t > 0, we have K^") (t) < W^"^ (t) , and for Tf,"\ < t < cr["^ , 
we further have from (4.18) that 

(4.28) K(-Ht) < W^-\t) = (W4")(r£\) - (t - rt],))^ . 

On the other hand, Proposition 4.3 and (4.17) with k replaced by A; — 1 
imply 

max \wt\s){- 00,0] A inf wi''\u)} 

= KW(ra) = w4")(aS-)V max^ ^ w(")(.)(-oo, 0] 



For t G [t^"\ , 0"^"^ ) , it follows from (4.18) and the above equality that 
K(")(t)= max|wj'^(s)(-oo,0] A inf W^|"^('u)| 

> max {wt\s){-oo,0] 

A inf VF^"^('u)A inf W^"'\u)} 

(4.29) = max {wt\s){-oo,0] A inf wi''\u)} 

A(M4")(rf_))-(t-ra))+ 
= A (W4"^(ra) - (t - ra)) + 

= (w4")(ra)-(t-ra))^ 
Equation (4.27) follows from (4.28) and (4.29). □ 



Remark 4.5. In light of (4.6) and Proposition 4.3, we have the charac- 
terization of r^"^ as 

^ = min{t > a'j^^\K^''\t) > W^''\t)} = mm{t > 4"V^"^(*) = 0}- 



(4.30) 

Because cr^-^ is the time of first arrival after t^"\ we in fact have 
(4.31) ?7(")(t)=0, rf)<t<4+\- 
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Evaluating (4.20) at ' and using Wg ' (a^ 
obtain 

(4.32) K(")(4"^) = VF^"^(4"^-). 

But (4.18) and Proposition 4.4 show that 



> K )(-oo,0], we 



(4.33) 
and so 
(4.34) 



in) 



<)(4")-), 



ak(")(4"^) = o. 



By contrast AK(")(t^"^) can be positive. Evaluating (4.20) at r^" and using 



(n) 



(4.17), we obtain 
(4.35) 

In conclusion, 



-in). 



(4.36) i^(")(i) = i^(")(4"V max >V^"^(s)(-oo,0], < t < r^"\ 



n). 



sea, 



(4.37) if(")(i) = (K(")(r(!\) - (t - r(::\))- 



4.2. Dynamics of the reference workload process. The evolutions of U^"'^ 
and are similar; the difference between them is asymptotically negli- 

gible. Before proving the properties ofU^^\ we provide a summary of these 
properties. The reader may work out the evolution of lA^'^^ and W*-"-* 

in Example 4.6 to follow along. This example appears in detail in [26]. 

Example 4.6. Consider a system realization in which 



Un 



Un 



(n) 
1 


= 1, 


(n) 
^1 


= 4, 




= 3, 




= 1, 


(n) 
2 


= 1, 


(n) 


= 4, 


rin) 
^2 


= 5, 


q{n) 


= 2, 


(n) 
3 


= 3, 


(n) 


= 2, 


tH 
^3 


= 1, 


qin) 


= 5, 


(n) 
4 


= 2, 


(n) 


= 1, 


rin) 


= 4, 


q[n) 


= 7, 


(n) 
5 


= 2, 


(n) 


= 1, 


rin) 
^5 


= 1, 


qin) 


= 9. 



«5 

Recall that 6a is a unit point mass at a. It is straightforward to compute 

ro, 0<t<l, 

{5-t)54-t, l<t<2, 

(5-t)54-j + 457-t, 2<t<5, 

(7-t)56-t + 4(57_t, 5<t<7, 

(ll-i)57-t + -5ii_f, 7<t<9, 

(^2(5_2 + 5i+52, t = 9, 
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0, 

(5-t)<54-t, 

(5 - t)6i.t + 4:67. 

{W-t)d-r.t, 
i8-t)6u-t, 
0, 

0, 

(5-t)54-t, 

(5-t)54-t + 4(57- 
(8-t)57-f, 

(9-t)57-f, 
(8-t)5ii_t, 

0, 
k 6i, 



0<t<l, 
l<t<2, 
2 < t < 4, 
4<t<5, 
5<t<6, 
6<t<7, 
7<t<8, 
8<t<9, 
t = 9, 
0<t<l, 
l<t<2, 
2 < t < 4, 
4<t<5, 
5<t <6, 
6 < t < 7, 
7<t <8, 
8<t <9, 
t = 9. 



Recall that K^"^^ is the amount of mass removed from the standard work- 
load W^l"^ to obt ain the reference workload [/("^ . To understand the process 
, we consider the dynamics of U^^\ In the absence of new arrivals, all 
atoms of ZY^") move left with unit speed. Moreover, the mass of the leftmost 
atom of U^""^ decreases with unit speed until it vanishes, corresponding to 
the work being done on the most urgent job in queue until it is served to 
completion [Proposition 4.8(i)]. However, if the leftmost atom of ZY^") hits 
zero, this atom is immediately removed from ZY^"-* [Proposition 4.8(ii), (v)]. 
This may be interpreted as reneging of a customer or deletion of a late cus- 
tomer from the system. When there is a new arrival at time t with lead time 
not smaller than the leftmost point of support of lA^^\t—), and this point of 
support is strictly positive, then a mass of the size v^2l^){t) lo'^^ted at ^^(i)^-^^ 

is added to U^"'\t—) [Proposition 4.8(iii)]. Similarly, if there is a new arrival 
and the leftmost point of the support of U^"^^ hits zero at the same time, 
then both of the above actions take place [(4.53) of Proposition 4.8(v)]. This 
is the case of a simultaneous new arrival and ejection of a late customer. The 
EDF system with reneging W^"^ shows the same behavior in all these cases. 
However, if a customer arrives to start a new busy period for U^^^ or, if at 
time t, there is a new arrival with lead time more urgent than the leftmost 
point of the support of U^'^\t—) (i.e., we have a "preemption"), then the 
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mass v^(n)(^fj associated with the new arrival is distributed in [Lj^(„j^^^,oo), 
or more precisely, on some atoms of W^"^ (t) located on this half-line, but it 

(n) 

is not necessarily located at the single atom L^(n)(^^y This possibility is de- 
scribed in Lemma 4.7 and Proposition 4.8(iv). In this respect, the evolution 
of U^^^ differs from that of W^"^ , for which all the new mass is always placed 
at the lead time of the arriving customer. Example 4.6 illustrates this. 

We now begin the rigorous study of As shown in Section 4.1, the 

time interval [0, oo) can be decomposed into a union of the disjoint intervals 
and > 0, such that J^(") = 1^^"^ - C/(") is non- 

increasing on (t^""*, o"^'^^-^] and nondecreasing on (<7^"\ ''"fc"'*]- In Lemma 4.7 
below, we analyze the behavior of ZY^"^ on the time intervals [t^"^,^^"^], 
k>l, while Proposition 4.8 describes the dynamics of U^^^ on the intervals 
(c7|."\ T^"^ ) , A; > 1. The section ends with Corollary 4.9, which describes the 
time evolution of the reference workload process 

We make use of the following elementary facts about the standard work- 
load. Since the interarrival times are strictly positive, A^(")(t) G {0, 1}, and 

(4.38) wPit) = wPit-) + AA(-\t)v^^l S^,., , t>0, 

^ ' A(")(t) 

which implies 

(4.39) A Wfit) = /\A^-){t)v^2lHt)^ * ^ 0' 

For any functions / and g on [0, oo) (taking finite or infinite values) such 
that whenever s <t and t — s is small enough, /(s) = f{t—) + t — s and 
g{s) = g{t—) + t — s, we have 

(4.40) ]imWP{s)[f{s),g{s)]=wP{t-)[f{t-),g{t-)]. 

This is true because the lead times of the customers present in the stan- 
dard system decrease with unit rate. Equation (4.40) remains valid if the 
closed intervals [/{■), g{-)] are replaced by either [/{■), gi-)), {f{-),g{-)] or 
(/(•), (7(-)). These facts will be used repeatedly in the following arguments, 
sometimes without mention. 



Lemma 4.7. Let k> 1. We have 



(4.41) C/W(t) = 0, rt\<t<a\:\ 

(4.43) U(-\aP){-oo,&l^^,„,^) = 0. 
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Proof. Equation (4.41) follows immediately from (4.6), (4.18) and Propo- 
sition 4.4. By (4.6), (4.34), (4.39) and the fact that A^(")(o-j,"^) = 1, we have 

AC/W(4")) = Al4")(af ) - Ai^W(4")) = 

and (4.42) follows. For y < , (4.5), (4.38), (4.34) and (4.33) imply 

ZYW(4"))(-oo,y] = K")(4"))(-oo,2/] -i^W(4"))]+ 

= K")(4")-)(-oo,y]-KW(4")-)]+ 

<[W^'^)(a(")-)-i^(")(4")-)]+ = 0, 
and so (4.43) also follows. □ 

(n) 

Lemma 4.7 shows that a\. begins a busy period for the reference system. 
Equation (4.30) imphes that [/(")(t) > for < * < , ^^'^ ^^us the 
intervals ) j k>l, are precisely the busy periods for the reference 

system. We analyze the behavior oiU^^^ during these busy periods. We start 
with the observation that, by (4.5) and Proposition 4.3, for t G (f t^"^ ) , 

(4.44) ^ 

WW(4")_)V max Wi")(.)(-oo,0])l . 

In what follows, given v £ M and any interval X C M, we will use 
to denote the measure in M that is zero on X'^ and coincides with u on 
T:u\j^{B) = v{BirM) for ah B£B{M). 

(n) (n) 

Proposition 4.8. For k>l and a), ' <t <tI , the following five prop- 
erties hold: 

(i) //A^(")(t) = and S(")(t-) >0, then 

(4.45) AJ^(")(t) = 0, 

(4.46) A[/(")(t) = 0. 

In this case, ifU^''\t-){E^''\t-)} > 0, then hothU'^'^\-){E^''\-)} andU^^'^t) 
decrease with unit rate in a neighborhood oft. On the other hand, if 

^W(t_){£;W(t_)} = 0, 
thenU(^\t)=wP{t)\^Ei-,^,^^^y 
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(ii) // AA(")(t) = and E^''\t-) = 0, then 

(4.47) (i-){0} = AK(") (t) = - a C/(") (t) 

and Z^W(t) = I (0,00). 

(iii) // A^(")(t) = 1 and lJ^j^^^ > E^'^\t-) > 0, i/ien ^^.^Jj /io/ds, 

A£;(")(t) >0 and 

(4.48) iYW(t)=wW(t-) + z;^t)m^L(") • 

(iv) //A^(")(t) = l, <-EW(t-), i/ien ^^.^5; /io/ds and 

(4.49) L%^^^<E(^^\t)<E(-\t-), 

(4.50) At/W(t) = z;^t)(i)> 

(4.51) Z^^"^(i)|(i^W(t_),oo)=^^"^(*-)|(£;(")(t-),oo), 

(4.52) ZY(")(t){^(")(t-)} >^Y(")(t-){^(")(t-)}. 

(v) // A^(")(t) = 1 and > E^'^Xt-) = 0, i/ien 

(4.53) u^n)(t)=U^-){t-) + v% 5^,^, _ZYW(t_){0}5o. 

^ A(")(t) 

Proof. Fix A; > 1 and t G (f^"\ t^""* ) • We start with the general obser- 
vation that, by (4.4) and (4.44), 

(t) = m.\ii{y\wP (t)(-oo, y] 



(4.54) 



>^5"^(4"^-)V max W^"^(s)(-oo,0]|, 



and because w|?^^(i) is purely atomic, the minimum on the right-hand side 

of (4.54) is obtained at the atom of {t) located at yo = E^'^\t). In 
particular, 

(4.55) >V^"^(t){s(")(i)} >0. 

We now consider each of the five different cases of the proposition. 

(i) Let a = £"(")(*-). By (4.4) and (4.5), for all s < t sufficiently near t, 

(4.56) wP{s){-oo,a/2]<K^''\s). 
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Also, for s £[t- a/2,t) sufficiently near t so that A(")(s) = holds 
[such s exist due to the assumption that AA^"'\t) = 0], we have 

(4.57) wP{t){-oo,0]<wP{s){-oo,a/2]. 

The last two relations show that W^"^(i)(— oo,0] < i^('")(t— ), and so, by 
Proposition 4.3, (4.45) holds. Equation (4.46) follows from (4.6), (4.45), 
(4.39) and the assumption A^(")(t) = 0. Because W^''\t) > [see (4.55)], 
decreases at unit rate in a neighborhood of t [see (2.6)-(2.8)]. In ad- 
dition, (4.6), (4.45) and the fact that, again by Proposition 4.3, X^"-* can- 
not increase on [c^"\ ''"^"^] except by a jump and hence is constant in a 
neighborhood of t, together imply that C/^"^ also decreases at unit rate in 
a neighborhood of t. Furthermore, the nature of the EDF discipline and 
(4.55) show that at t, the standard system is serving a customer with lead 
time no greater than E^^\t). Combining the above properties with the 

fact that ^^"■'(*)|(£;(")(i),oo) = ^s"''(*)|(£(")(t),oo) (4-9), we conclude that 
if Z^(")(i-){S(")(t-)} > 0, then U^"'\-){E^''\-)} decreases with unit rate in 
a neighborhood of t. On the other hand, if U^^''\t—){E^'^\t—)} = 0, then 
since jumps up at t. Indeed, in this case, 

>v[r)(t)(-oo,s(")(t-)] = wJ'^(t-)(-oo,E(")(i-)] = K(")(t-) = irW(t). 

This means that 

E(")(t) = min{y G R\W^^\t){-oo,y] > ^("^(t)} 

= min{y > sW(t-)|wJ*^(t){y} > 0}. 

It follows that 

(4.58) Wj'^(t)(^(")(t-),^(")(t)) = 0. 

Using the definition of S(")(t), (4.46), (4.9), the assumption Z^(")(t-){£;(")(t-)} 
= 0, the assumption AA(")(t) = 0, and (4.58), we obtain 

ZY(")(t)[E(")(t),oo) = /7(")(t) = [/(")(*-) 

= ZY(")(t-)[E(")(t-),oo) 
= ZY(")(t-)(^(")(i-),oo) 

= wP{t){E^'^\t-),^) 
= wP{t)[E(^Ht),^). 



28 KRUK, LEHOCZKY, RAMANAN AND SHREVE 

From (4.9) we see now that ^^"HO = W?^(0|[£;W{t),oo)• 
(ii) By (4.30), (4.4) and (4.5), for s £ (4"^*) we have 

(4.59) wJ'^(s)(-oo,s(")(s)] >i^(")(s). 

As st* in (4.59), by (4.40), (4.38), and the case (ii) assumptions AA^'^\t) = 
and £;(")(*-) =0, we get 

(4.60) Wlr^(t)(-oo,0] = wJ"^(t-)(-oo,0] >K(")(t-). 
When combined with Proposition 4.3, this imphes 

(4.61) K'^''\t) = K^''\t-)VWP{t){-oo,0]=wP{t){-oo,0]. 
By (4.4) and (4.5), for s G (4"\i)> 

Kin) (^){^(n) (^)l ^ W (^s){-00, ^("^ (s)] - K^") (s). 

Letting s'[t, invoking (4.40), (4.60) and (4.61), and recalhng the assumption 
£:W(i-) =0, we obtain 

(4.62) ZY(")(i-){0} = wP{t-){-oo,0] - K^'^Xt-) = K^''\t) - K^^'^t-), 

and the first equahty in (4.47) follows. The second equahty in (4.47) follows 
from (4.6), (4.39) and the assumption AA^^\t) = 0. Moreover, by (4.5) and 
(4.61), for every y G M, 

(4.63) 

= W?^(0|(o,oo)(-oo,2/]. 

(iii) Let a = E^"-\t-). We can deduce (4.45) from (4.56) and (4.57) as 
in (i), with the only difference that now (4.57), for s < t sufficiently close 
to t such that A^"-\t-) = A(")(s), follows from the fact that L^^l)^^-^ > 0, 
since this implies that the work for the system associated with the customer 
arriving to the system at time t does not contribute to VV^"^(t)(— oo, 0]. Next, 
let y < a, let e = {a — y) /2 and note that by assumption, L^^(„) > a > y + e. 
Thus, for s < t, s sufficiently close to t [so as to ensure that A(")(t-) = 
AW(s)], we have wP{t){-oo,y] < {s){-oo,y + e] < iv:W(s), where the 
last inequality uses (4.5) and the fact that y + e< E^'^\t-). Letting s ft, 
we obtain W^"^(t)(— oo,y] <K^^\t—), which, together with (4.45), shows 
that y < E^''\t). Thus, E(")(t-) < £;(")(*) or, equivalently, A£;(")(i) > 0. 

We now turn to the proof of (4.48). Equation (4.5) implies 

(4.64) ZY(-)(t-)(-oo,y] = [Y^f {t-){-^,y] - K^^\t-)]^ . 
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Indeed, for any y such that VV^^"^ (t— ){y} = 0, (4.64) follows from (4.5), in 
which t is replaced by s < t, by taking s '[ t. However, the family of sets 
{—oo,y] with y^g^\t—){y} = forms a separating class in ;S(]R), and so 
(4.64) holds for all y. Moreover, using (4.38), (4.45) and (4.5), we see that 

ZY(")(t)(-oo,2/] 

(4.65) 

= [WP{t-){-^,y]-K(^\t-)+v% 6^,., (_oo,y]] + . 

When combined with (4.64), this shows that 

(4.66) ^(")(t)(-oo,y] =^W(t-)(-oo,y], y < L^^l ^^^{t) . 

On the other hand, if y > ^J(l)(j)(i)> then y > ^(")(t-) and (4.64) becomes 

^(")(t_)(_oo,y] = W^")(t-)(-oo,y] -i^(")(t-). 
From (4.65), we now have 

(4.67) ^W(t)(-oo,2/] =U(-\t-){-oo,y] + v%^^^, y > L%^^^. 

When combined, (4.66) and (4.67) prove (4.48). 

(iv) We have -^^(i)^^) > 0, and so (4.45) holds by the same argument as 

in case (iii), but now with a = -^^(Ij^^y The assumptions -^^^{l)^^^ < E^^\t—) 

and AAW(t) = 1, along with the relations (4.38), (4.5), (4.45) and the def- 
inition of E^'^\ imply that 



>wf\t-){-oo,E^''\t-)] 

>K(")(t-) 

= K(")(t). 

Invoking (4.5) again, this shows that U^^\t){—oo,E^^\t—y\ > 0, which im- 
plies E('^)(t) < E(")(t-). Now, let y <a = l'-^I^^^^ and let e = {a - y)/2. 

Then, combining (4.38), the inequalities y + £ <a< E^'^\t—) and (4.45), we 
obtain 

wP{t){-oo,y]<W^^\t-){-oo,y + e]<K(^Ht-)=K(^\t). 

This shows that y < E^'^\t), which proves (4.49). In addition, by (4.6) and 
(4.45), we have 

C/(")(t) = VF(")(t)-/^W(t) 
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(n) 



and (4.50) follows. Furthermore, since £'(")(*) < E^"-\t-) by (4.49), the re- 
lations (4.9), (4.38) and the assumption ^^(i)^^) < E^^^^t-) imply 

= ZY(")(t-)|(^(„)(i_)^^). 

This estabhshes (4.51). 

Finally, to prove (4.52), we will consider two cases. 
Case I. £;(")(t) <S(")(t-). 
By (4.9), we know that 

^YW(t){£;W(t-)} = H;W(t){^(«)(t_)}. 

In turn, when combined with (4.64) and the definition of E^'^\ this shows 
that 

i^W(t_){£;W(t_)} 

= (t-) (-00, (t-)] - (t-) (-00, (t-)) 
= (t-) (-00, (t-)] - (t-) 

- K"Ht-)(-oo,i5;(")(t-)) -i^(")(t-)]+ 
<wP{t-){-oo,E^^\t-)]-WP{t-){-oo,E^^\t-)) 
= wP{t){E(^\t-)} 
= ZY(")(t){^(")(t-)}, 

and so (4.52) holds. 

Case II. #")(t) = S('^)(t-). 

By (4.5), (4.38), (4.45), (4.64) and the definition of ^("), 
UH (t)} = (t)(-oo, (t)] 

= Wi"^ {t){-oo, (t)] - (t) 

= wf (t-)(-oo,i?(")(t-)]+^;^t)w-^^"^(*-) 

(n) 



:wW(t-)(-oo,£;W(t-)] + ^5;i,^^) 
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(n) 

A(")(t)' 



which estabUshes (4.52) in this well. Since S(")(t) <E('^\t-), the 

two cases above are exhaustive, and so (4.52) is proved. 

(v) Equation (4.63) holds by the same argument as in (ii), but where now 

the equality in (4.60) follows from the fact that -^^^(l)^^-) > 0. Let u[^\t) = 
ZY(")(t-) + 7;^"!) (5^(„) -ZY(")(t-){0}5o. We want to show that ZY(")(t) = 

U[''\t). By (4.10), ZY(")(t) andZY(")(t-) are supported on (0,oo) and [0,oo), 
respectively. Thus, 

(4.68) ^/(")(t)(-oo,y]=ZYj")(t)(-oo,y]=0, y < 0. 

By (4.9) and the fact that £;(")(*-) = 0, we have ZY(")(t-)|(o,oo) = >V5"^(t-)|(o,oo) • 
The last two statements, along with (4.38), (4.63) and another application 
of (4.9), show that 

Z^f )(t)|(0,oo) =W^"^(t-)|(0,oo) 

^ ' A(")(t) 

= W^")(t)|(0,oo) 
= ZY(")(t)|(0,oo). 

This, together with (4.68), shows that ZY(")(i) = ZY|"^(i). □ 

The last result of this section concerns the evolution of Despite 
the different ways in which arriving mass is distributed in the system with 
reneging and the reference system, in both systems one can keep track of the 
total mass in system by beginning with the arrived mass (which is the same 
in both systems), subtracting the reduction in mass due to service (which 
occurs continuously at unit rate per unit time whenever mass is present), and 
subtracting the mass that has become late and been deleted. In particular, 
a simple mass balance shows that 

(4.69) W'-^Ht) = y(")(^W(t)) - ^*I|VK("){.)>0} ds - <Ht), 

(n) 

where we recall that R^y is the total amount of reneged work in the reneging 
system, which admits the representation (2.18), R^\t) = X]o<s<i 
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for all t G [0,00). We now show that the following analogous relation holds 
for the reference workload: 

(4.70) C/W(t) = - ^*l{f;(.,(.)>o} ds - 4'\t), 
where 

(4.71) <)(t)^ Yl ^^"^(^-){0}. 

0<s<t 

Also, for notational convenience, we set R^\o—) = R^'\o—) = 0. 

Corollary 4.9 
and hence 

(4.72) 

where, for t>0, 
(4.73) 

Proof. For t>0, let U^"-\t) be equal to the right-hand side of (4.70). 
By (4.41) of Lemma 4.7, we have U^"-\0) = = U'^"-\0). Moreover, for every 
A; > 1, by Lemma 4.7 and the definition of o'^"'*, it follows that U^^\t—) = 
;7W(t) =0 and AT/W(A(")(t)) = for i G (t^"\,4"^), C/W(4"^-) = and 
AC/(")(ct!."^) = AT/(")(^(")(f7^)). When compared with the right-hand side 
of (4.70), this shows that and are both flat on (r^"\ , cr|."^ ) , with an 
upward jump at cr^"^ of size AfW(A(")(4"^)). Thus, to prove the corollary, 
it suffices to show that the increments of C/^"^ and C/^"-* on the intervals 
(<t["\t^"^], A: > 1, coincide. 

Fix /c > 1. We first show that 

(4.74) A[7(")(Tf ^) = A;7(")(Tf ^). 

Equality (4.30) shows that there cannot be an arrival at time for such 

(n) 

an arrival would have a positive lead time and hence increase Wf' with- 
out increasing (see Proposition 4.3). In other words, /S.A^'^\t'^'^) = 0. 
Because there is no arrival at r^"^ , the measure- valued process W^"^ is con- 
tinuous at T^'^ . Taking the limit in (4.5) as t'\T^\ we obtain 

Z^^"^(rf^-)(-oo,0] = [w(")(rf))(-oo,0]-irW(rf)-)]+ = Ai^W(r(")), 



For every t>0, equation (4- 70) holds. Moreover, 

ijj^ (t) = J I{C7(") (s)=0} ds. 
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where the last equahty is a consequence of (4.36). However, (4.10) impHes 
that ZY(")(Ti"^-)(-oo,0) < hm^^ („) ^/(")(t)(-oo,0) = 0, so AC/(")(t^"^) = 

_^(n)(^C«)_^)|Q| ^ _AiC(")(r^")). From (4.6) and the continuity of 

at r^"\ we see that -AirW(T^"^) is also equal to AC/(")(t^"^), and (4.74) 
is proved. 

We next show that Ac7(")(t) = AC/W(t) for t G (4"^ ^i"'')- If E^'^Ht-) > 
0, then the definitions of E^'^') and f/'^"-', and statements (i), (iii) and (iv) of 
Proposition 4.8 show that 

AC/(")(t) = AC/(")(t) = AF(")(A(")(t)). 

On the other hand, if E^^\t—) = 0, then properties (ii) and (v) of Proposi- 
tion 4.8 and the definition of f/^"^ show that 

AC/W(t) = AC7(")(t) = AA(-\t)v%^^^-U'^-\t-m. 

Now, let 5^") be the (random) set of times s > for which U^'^\s) > 
and at least one of the following three properties holds: 

AA(")(s) >0, 

£;W(s-) = 

or 

^yW(s-){£;W(s-)} = 0. 

Suppose ?7(")(s) > 0. If ^(")(s-) = 0, then the fact that ^(")(s) > by 
(4.10) implies A£;(")(s) > 0, while if ^YW(s-){^(")(s-)} =0, the defini- 
tion of £;(")(s) implies that A(ZY(")(s){^(")(s)}) > 0. Thus, the set 5^") 

is countable, and on the set {s G {al^\t) :C/(")(s) > 0} \ S^"\ the process 
[/^"^ decreases with unit rate by Proposition 4.8(i). Therefore, the total 
amount of this decrease on any time interval of the form {a^\t) equals 
/^(") I{c/(")(s)>o} which coincides with the absolutely continuous part of 

C/('^)(t) - U^^'HaP-) on the same interval. This concludes the proof of 
(4.70). 

Adding and subtracting t to (4.70), by the definition (2.6) of the netput 
process iV(") and the nonnegativity of we obtain 

(4.75) C/W(t) = ivW(t)+ ri^^(„,(^)^,^d.-i?W(t), 



while substituting (4.15) and (2.8) into (4.6), we have 

(t) = (t) + 4"^ (t) + xi"^ (t) - K^^^ (t) 
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for t > 0. On the other hand, we know that 

/ \u(")is)>o}dIs'\s) = and / I|^{„)(^)>o} ^^^-"^(•s) = 0, 

^[0,oo) J[0,oo) 

where the former equahty holds because W^'^ > C/^"^ by (4.8), and 

increases only at times when VFg"^ is zero, while the latter holds by (4.16). 
From the last three displays, we conclude that 

(4.76) / I^^(„)(,)>o}d4")(s)= / I^^(„,(,)>ojdKf (s). 

J[0,oo) J[0,oo) 

On the other hand, since U^"'\s) = implies AA^'^\s) = 0, from properties 
(i) and (ii) of Proposition 4.8 and the fact that is a pure jump process 
with AR'i^\t) =iY(")(t-){0}, it follows that 

J[0,oo) J[0,oo) 

(n) (n) 

Together, the last two equalities imply Rjj = , which, when substituted 
into (4.70), yields (4.72). □ 

5. The reneging system. In this section we bound the difference in work- 
load between the pre-limit reference and reneging systems — Lemma 5.2 pro- 
vides a lower bound, while Lemma 5.6 provides an upper bound. The proof 
of the upper bound uses an optimality property of EDF that may be of 
independent interest. 

Theorem 5.1. Let ir be a service policy for a single-station, single- 
customer- class queueing system with reneging such that the customer arrival 
times to this system do not have a finite accumulation point. Let R-jrit) be 
the amount of work removed from this system up to time t due to lateness. 
Let Rwit) be the amount of work removed due to lateness up to time t from 
the EDF system with reneging and the same interarrival times, service times 
and lead times as in the former system. Then for every t>0, we have 

(5.1) Rwit)<Rn{t). 

The proof of Theorem 5.1 is deferred to the Appendix. The related fact 
that the EDF protocol minimizes the number of late customers in the G/M/c 
queue was proved in [29] , and the main idea of our proof is similar to that 
of [29]. However, our argument is pathwise and the only assumption on 
the distribution of the system stochastic primitives that we impose is that 
customer arrivals do not have a finite accumulation point. This assumption 
is clearly satisfied almost surely by a GLjGjl queue. 
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5.1. Comparison results. In this section, we establish bounds on the dif- 
ference between the processes C/^"^ and W^"'h In Section 6.1, this difference 
will be shown to be negligible in the heavy traffic limit. We start with Lemma 
5.2 showing that W^'^'^ < U^"'\ which implies that R^^^ < Ry^^ (see CoroUary 
5.3). 

In the proofs of these results, we will make frequent use of the observation 
that, by (4.69) and (4.70), 

H^W(t)_[/H(t)= [\ruinu,.^o}ds- f'hw('')(s)>o}ds 



(5.2) 

+ R^-\t) - R^-\t) 

for t G [0,oo). 

Lemma 5.2. For every t>0, we have 
(5.3) W^(")(t) (t). 



Proof. Let 

(5.4) r = min{t >0:W^''\t) > U^''Ht)}. 

If r = +00, then (5.3) holds. Assume r < +oo. In this case, we claim that 
the minimum on the right-hand side of (5.4) is attained. Indeed, (5.2) and 

the fact that and R^^ are pure jump processes show that the only way 
that W^"^^ — [/'•"^ can become strictly positive is via a jump. Thus W^"'\t) > 
[/^(t). Since >V(r)(-cx),0] = ^/(t)(-oo, 0] = (in fact, this equality holds 
for any time t), this means there must exist a y > such that 

(5.5) W(")(r)(y,oo)>^(")(r)(?/,oo). 
Let 

(5.6) ro = inf{te [0,r]:>v(")(t)(y + T-t,oo) >ZY('')(t)(?/ + r-t,oo)}. 

By (5.5), the above infimum is over a nonempty set. Lemma 4.7 and Propo- 
sition 4.8 imply that the only difference in the dynamics of W^"^ and U^'^^ 

(n) (n) 

is that the arriving mass is concentrated at in the case of the EDF 

system with reneging and distributed in [L^^^oo) in the reference system. 
On the other hand, in both systems at time t € [0, r] , no mass leaves the 
interval (y + t — t,oo) due to lateness. This implies that W^"'\t){y + r — 
t,oo)-U^'^\t){y + r — t, oo), t G [0, r], has no positive jumps and therefore 



(5.7) W(")(t)(y + T - To,oo) =Z^(")(t)(y + r - ro,oo). 
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By (5.5) and (5.7), tq < r. Thus, there exists t G (ro,r), where t — tq is 
arbitrarily small and 

(5.8) W^"'\t){y + T-t,oo)>U^''\t){y + T-t,oo). 

However, we claim that (5.7) and (5.8) imply that for all t G (ro,r), where 
t — Tq is small enough, it must be that 

(5.9) W(")(0(0,y + r-t] >0, 

(5.10) Z^(")(i)(0,?/ + r-t] =0. 

Indeed, if (5.9) is false, then the left-hand side of (5.8) is equal to W^^\t), 
and consequently decreases with unit speed as long as it is nonzero in some 
time interval beginning with tq. Similarly, if (5.10) is false, the right-hand 
side of (5.8) is constant on some interval beginning with tq. In both cases, 
due to (5.7), (5.8) cannot hold for t G (to,t) with t — tq arbitrarily small. 
But (5.8)-(5.10) yield W'^'^\t) > U^'^\t) for some t < r, which contradicts 
(5.4). □ 



Corollary 5.3. For every t>0, 

(5.11) <)(t)<<)(t). 

(n) 

Moreover, for k>l and t>ai, , 

(5.12) RP{t) - <Vf -) < R%\t) - <Vf - 



Proof. Lemma 5.2 and (5.2) imply that for < s < t, 

(5.13) (4") (t) -<)(.)) - {R^^\t) - R%\s)) < [/W(.) -T^W(.). 

Substituting s = into (5.13) and using the fact that r'^\o) = R^^\o) = 
;7(")(0) = VF(")(0) = 0, we obtain (5.11). Likewise, for < s < < t, tak- 
ing limits as s tends to cj^"^ in (5.13), and using the fact that U^"'\a^^ —) = 
W^'^\a^^ —) = 0, which follows from (4.41), Lemma 5.2 and the nonnega- 
tivity of W'^"), we obtain (5.12). □ 

The proofs of Lemma 5.2 and Corollary 5.3 show the following more 
general (and intuitively obvious) fact: if all customers in the EDF system 
with reneging get larger deadlines, this results in a larger workload at every 
time t and a smaller total amount of mass removed from the system due to 
lateness in the time interval [0, t] . 

We now establish an inequality between the frontiers in both systems. 
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Lemma 5.4. For every t>0 such that U^'^\t) > 0, we have 
(5.14) E(")(t)<F(")(t). 

Proof. Subtracting (4.5) from (4.6), we see that for any y £ M, 

ZY(")(t)(y,oo) = l^(")(t) - KW(t) - [wPit)i-oo,y] - i^(")(t)]+ 

(5.15) 

<>v(")(t)(y,oo). 

Now, assume that for some t we have F^"'\t) < E^'^\t). In this case, 
ly (t) > W^") (t){C7W (t)} + (0 (F(") {t),oo) 
= (t){C(") (t)} + V(") (t) (f(") (t), oo) 

> W(")(t){C(")(t)} + v(")(t)[^(")(t),oo) 

(5.16) 

> W('')(t){C('^)(t)} + wf'^(t)[^(")(t),oo) 

> w(")(t){C(")(t)} +ZY(")(t)[^(")(t),oo) 

>c/^"Ht), 

where the second hne follows from the fact that none of the customers in 
the EDF system with reneging that have lead times at time t greater than 
F^^\t) has received any service up to time t, the second-last inequality 
follows from (5.15) and the last line holds due to the equality U^'^\t) = 
Z^('^)(i)[^('')(t),oo). When combined with the assumption that C/(")(t) > 0, 
this implies that W^'^\t) > 0. This, in turn, implies that >V(")(t){C(")(t)} > 
because the residual service time of the currently served customer is strictly 
positive. Thus, the last inequality in (5.16) is strict, which contradicts (5.3). 
□ 

Let D^"'\t) be the amount of work deleted by the EDF system with 
reneging in the time interval [0,t] that is associated with customers whose 
lead times upon arrival were smaller than the value of the frontier at the 
time of their arrival. In the proof of the next lemma, we will make use of 
the elementary fact that by the definition of F^") we have 

(5.17) F(")(ti) - (t2 - ti) < F(")(t2), <ti< t2. 

Lemma 5.5. For every t>0, 

(5.18) C/(")(t)-VI/(")(i)<Z?(")(i). 
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Proof. If t G [r^"\,4"0 for some /fc > 1, then (t) = by (4.41). 
Thus, by (4.19), it suffices to prove (5.18) on [f^""* , t^""* ) for every A; > 1. Let 
A; > 1. Suppose that (5.18) is false for some t G ''"fc"'* ) • Let 

(5.19) T = min{t G [^7^"^-^^" W " W^'^Xt) > Z)(")(t)}. 

We first argue that the minimum on the right-hand side of (5.19) is attained. 
Indeed, by (5.2) and Lemma 5.2, it is clear that C/^"^ — Vl^^") cannot increase 
except by a jump that is due to lateness in the EDF system with reneging. 
Thus, we have >V(")(t-){0} > and 

(5.20) C/(")(r)-W^(")(^)>^^''H^)- 

Also, (4.41), (4.42) and Lemma 5.2 imply that [/("^(cj^''^) = AC/(")(4"^) = 
AVF(")(4"^) = VF(")(4"^), so < -r- In particular, (5.19) implies 

(5.21) ?7(")(t-) - M^(")(t-) < D(")(t-). 

(n) (n) 

Let ko be the index of the customer arriving at time <t^ , that is, Sj^^ = 

(n) 

o"^ . Let ki > ko be the index of a customer who reneges in the reneging 
system at time r. There must be such a customer, and there may in fact be 
more than one such customer. The amount of work associated with all such 
customers at time r is W^"'\t—){0} , and we seek to show that this work is 
bounded above by AD(") (r) . We have sj^^ G [^"^ , r) and lJ^"^ - (r - sj^^ ) = 
0. The subsequent analysis is divided into two cases. 

Case I. For every customer ki chosen as just described, assume there is a 
customer i arriving in the time interval [c^"^ , S^"^] who is at least as urgent as 
customer ki when customer ki arrives but whose associated mass in the ref- 
erence system is at least partly assigned so that upon the arrival of customer 
ki , this mass is to the right of L^""* . In other words, £ £ [ko, ki] , L^""* — (-S^."^ — 
5f )) < 4"^ and AWW(sf )){lV^} > AZYW(5f ))[l{"\4';) + 5^';) - 4"^]. 

In this case, A^(")(sf ^)(L^"^ + sj^'^ - \oo) > 0. Indeed, by Lemma 4.7 
and Proposition 4.8 (iv) (describing the only case in which part of the mass 
of a new customer is distributed by the reference workload to a point other 
than its lead time) A^7(4"^) = v'f^ and AW(")(5f ^)(-oo, 4"^) = [see 
(4.42), (4.43), (4.49), (4.50) and (4.4)]. Let s > 4"^ + 4f " ^e"^ satisfy 

AZY^"'^(5'^"^){s} > 0. Such a point s exists since the measure U^^\S^"''^) is 
discrete. 

If ^ > fco (e.g., £ = ki), then, by (4.51) in Proposition 4.8(iv) and Lemma 
5.4, we have s < ^("^(sj"^-) < F(")(5^"^-) < F(")(4"^). Thus, by (5.17), 
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U£ = ko, then, because ZY(")(5£^){s} > 0, we have wj'^(5^"^){s} > by 
the definition of In this case W^") (5i"^){s} = 0, because = on 
[tS,^^"^) by (4.41) and Lemma 5.2, so = = 4o^<5^(n) 

feQ 

and s > + sj^^ — S^^^ > L^j^ by the definitions of £ and s. Thus, a 

(n) 

customer with lead time equal to s at time Sj^^ has already been in service 

in the EDF system with reneging, so ^i"^ + 5^"^ - < s < and 

consequently, by (5.17), L^^'^ < FW(5g)) - (k:^ - 5^^) < ^^(^fV 

Thus, regardless of the value of £, L^^^ < F^^\S^'^). In other words, un- 
der the Case I assumption, every customer ki who becomes late at time r 
in the EDF system with reneging arrived with initial lead time smaller than 
the value of F^"^^ at the time of its arrival. The work associated with these 
customers deleted at time r is AD(")(r). We conclude that W(")(t-){0} = 
AD(")(r). However, by (5.2), we have A(?7(") - W^("))(r) < W(")(r-){0}, 
and so A(C/(") - 1^("))(t) < AZ)(")(t). This, together with (5.21), contra- 
dicts (5.20). 

Case II. For a customer ki chosen as described above, assume that every 
customer £ arriving in the time interval [f^,""*, 5*^"^] who is as least as urgent 
as customer ki when customer ki arrives has all its associated mass ini- 
tially assigned in the reference system to the interval (0,L^"^ -|- — S^^^] 
upon arrival. Customers £ who are less urgent then ki must have lead times 
satisfying L^"^ > L^"-* -|- S"^"^ — S^"'\ and hence the mass brought by such 

customers must be initially assigned to the half-line {L^^^ + — S^"'\oo) 

in both systems. Then for every t G Wk^\ S^^], we have 

(5.22) W(")(t)(0,4;) - (t - 5^;^))] < ^(")(t)(0,4;) -{t- 5^;^))], 

as we now explain. Under the Case II assumption the arrival of new mass 
is the same on both sides of (5.22). Furthermore, disregarding lateness and 
new arrivals, both sides of (5.22) decrease at unit rate so long as they are 
nonzero. Finally, by (5.12) the amount of late work removed from the EDF 
system with reneging in the time interval , t] is greater than or equal to 
the amount of late work removed from ZY^") in this time interval. Therefore, 
(5.22) holds for every t G [(tI!'\ sj^^]. 

We claim that (5.22) in fact holds for all t G [<t^"\''")- Suppose this is not 
the case. Let 

r?^inf{t G [5i:\r)|w(")(t)(0,4':) - (t-Sj^^)] 

(5.23) 

>zY(")(t)(o,4';)-(t-5(f)]}. 
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The strict inequality in (5.23) can occur only because of an arrival at time t 
which brines mass to the interval (0, lJ."^ -{t- S"^"^)] under the W^") mea- 
sure but not under the U^"'^ measure. The arrival at time ki does not have 
this property because the Case II assumption applies to i = ki. Therefore, 

Also, for t G [5^"\t), 

(5.24) >vW(i){4n)_(i_5W)}>o, 

because the customer ki is present in the EDF system with reneging at time 
t. By (4.4), (5.24) and the definition of r?, we have ^(")(t) < L^"^ -{t- 5"^"^) 

for t e r?). Thus, £;W(t-) < L^^^ -{t- S^"^) for t G (S^"\r/]. We argue 
that this implies that the amounts of mass arriving in both the EDF system 
with reneging and the reference workload at any time t € {S^^ ,7]] with lead 

times upon arrival less than or equal to L^"^ — (^ — ^hi^) the same. Indeed, 
Proposition 4.8, especially (4.51), implies that no mass arriving at time t 
with lead time smaller than E^'^\t—) in the EDF system with reneging is 
distributed to lead times greater than E^'^\t—) by the reference workload. 
Also, Proposition 4.8 (iii) and (v) imply that the mass arriving at time t with 
lead time greater than or equal to E^'^\t—) is distributed in the same way 
by the EDF system with reneging and the reference system. By the same 
argument as in the case of t G ['^^if' : ^t^^]^ conclude that (5.22) holds for 
t G [S'|,"\r/], which contradicts the definition of r]. We have shown that (5.22) 
holds for t G [4"\r). 

Letting t\ t in (5.22) and using the fact that L^"^ - (r - 5^"^) = 0, 

we get W(")(r-){0} < ZY(")(r-){0}. Thus, by (5.2), A(?7(") - VF("))(r) = 
W('")(r-){0} -^(")(r-){0} < which, together with (5.21) and the fact 
that is nondecreasing, contradicts (5.20). □ 

For the sake of the next proof, we define a sequence of auxiliary hybrid sys- 
tems (with the same stochastic primitives as in the case of the EDF systems 
described in Section 2.2) as follows. The hybrid system gives priority to the 
jobs whose lead times upon arrival are smaller than the current frontier 
in the corresponding EDF system with reneging. In other words, for each k, 
the kth customer arriving at the hybrid system joins the high-priority class 
if and only if 

(5.25) 4"^ <F(")(4")). 

The system processes high-priority customers according to the FIFO service 
discipline. When the priority class empties, the system goes idle until either 
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another high-priority customer arrives and the system resumes service in the 
manner described above, or the corresponding EDF system with reneging 
finishes serving the customers who have received priority in the hybrid sys- 
tem. Here, we are using the fact that the high-priority customers leave the 
hybrid system before they leave the EDF system with reneging, which is a 
consequence of the optimality of the EDF discipline established in Theorem 
5.1. (We have slightly abused the terminology here, identifying the kth cus- 
tomer in the hybrid system with the corresponding customer from the EDF 

(n) (n) 

system with reneging, while, formally, only the random variables uj^ , vj^ 

and L^""* associated with these customers are the same.) Whenever the EDF 
system with reneging finishes serving a batch of customers who have received 
high priority in the hybrid system, both systems then serve the low-priority 
class using the EDF discipline until the next high-priority customer arrives. 
In both systems, if a customer is present when his deadline passes, he leaves 
the queue immediately, regardless of his class. The measure- valued workload 

(n) 

process associated with the hybrid system will be denoted by Wjj . 
Lemma 5.6. For every t>0, we have 

yl(")(t) 

(5.26) < Y: -i"^(wW(5i")-)(0,FW(5f)) + 4")-4"))+ 

k=l 

^\4"'<F(")(sf^)}- 

Proof. By Lemma 5.5, it suffices to show that D^"'\t) is not greater 
than the right-hand side of (5.26). By Theorem 5.1, D^'^\t), the amount 
of unfinished work associated with customers who arrived with lead times 
smaller than F^'"'^ and were deleted in the time interval [0,t] by the EDF 
system with reneging, is not greater than the unfinished work associated with 
these customers and deleted by the corresponding hybrid system. Note that 
the customers with lead times satisfying (5.25) form a priority class in both 
the EDF system with reneging and the hybrid system, and so their service 
is not affected by the presence of other customers. Furthermore, unfinished 
work associated with deleted customers who arrived with lead times greater 
than or equal to F^"'^ is the same in both systems. 

For each k, if (5.25) holds, then the kth customer of the hybrid sys- 
tem belongs to the high-priority class. Moreover, if, for some i<k, l'i"^ < 
fW(5{"^), then, by (5.17), l}"^ - (5^"^ - ^f"^), the lead time of the Ith 
customer at time S'^\ is smaller than F^"'\S^'^). Thus, if (5.25) holds, the 
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kth. customer waits at most W^/^\s^^ — ){0, F^^\S^j^^ )) time miits before he 
starts receiving service. (His waiting time may actually be smaller because 
some of the high-priority customers in queue who have arrived before him 
may renege before they are served to completion.) We have 

(5.27) w{;'^(5^")-)(0,F(")(5^"))) < W(")(sf^-)(0,F(")(5^"))), 

because, in both systems under consideration, the arrivals with lead times 
smaller than F^""^ and the corresponding work associated with them are the 
same, the server serves these customers with rate 1 as long as they are present 
in the system, but, by Theorem 5.1, the amount of unfinished work associ- 
ated with these customers and deleted by the EDF system with reneging is 
not greater than the work deleted by the hybrid system. Thus, if (5.25) holds, 
the time required for the hybrid system to fully serve the fcth customer is at 

most W(")(5^"^-)(0,F(")(5^"^)) + 4"^ Therefore, under assumption (5.25), 
the unfinished work deleted by the hybrid system due to lateness of the kth 
customer is at most 4"^ A (W(")(5^"^-)(0, F(")(5^"^)) + 4"^ - Thus, 
the amount of work associated with high-priority customers deleted by the 
hybrid system up to time t is bounded above by the right-hand side of (5.26). 
□ 



6. Heavy traffic analysis. In Sections 6.1 and 6.2, respectively, we iden- 
tify the heavy traffic limit of the scaled workload and the scaled reneged 
work in the reneging system. In both cases, this is done by first considering 
the reference system, which is easier to analyze, and then using the bounds 
derived in Section 5.1 to show that the limits in both systems coincide. For 
the heavy traffic analysis of the reference system, we will find it useful to 
introduce the following scaled quantities: 

Jn 

(6.1) 



and, for every Borel set S C M, 

(6.2) Z?(")(t)(S) = ^ZY(")(nt)(V^5). 



Jn 

Also, define 

(6.3) U* = ^{Ws) and U* {■) =U* = ^{Ws){n). 

6.1. Proofs of main results concerning the workload. 



6.1.1. Proof of Theorem 3.4- In Lemma 6.1, we use the continuity prop- 
erty of the mapping <I> established in Lemma 4.1, along with the character- 
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ization of the heavy traffic hmit of the workload measure- valued process in 
the standard system, to identify the heavy traffic limit of the workload in 
the reference system. Let Aj|^(o) : -D[0, oo) — )• D[0, oo) be the mapping defined, 
for every (p £ D[0, oo) and t > 0, by 



(6.4) AH(o)(</')(t) = sup {ct){s) - H{0))+ A inf 0(n 



se[o,t] 



If (p is nonnegative, then by Theorem 1.4 from [25], Aj:/(o) ((/>) is the function 
in D[0, oo) obtained by double reflection of (/> at and H{0). In other words, 
^H{o) {4') takes values in [0, Hifi)] and has the unique decomposition 

(6.5) A^(o) ((/>) = - K+ + K„ , 

where k± are nondecreasing RCLL functions satisfying k±(0— ) = and 

{kH(Q)ms)<Hmdi^+{s) = ^, 

0,00 

(6.6) 

I{A^(O)(0)(5)>O} dK-(s) =0. 

[0,oo) 



Lemma 6.1. The process U* satisfies 
(6.7) U*=AHio)iW*s) 

and has the same distribution as W* . Moreover, U^"^^ =^hl* = $(W^) and 
fj{n) ^^r* asn^oo. 

Proof. By the definition of U* and $ given in (6.3) and (4.2), respec- 
tively. 



U*{t) = '^{W*s)iR)it) = Wsit)- sup Ws(-oo,0]A inf VF^(n 



s6[0,i] 



ue[s,t] 



t > 0. 



Since (3.1)-(3.3) imply W^{t){-oo,0] = {W^ (t) - H {0))+ for every t > 0, this 
shows that U* = Ah(^o)0^s)- ^^'^ characterization of Wg given at the end 
of Section 2.4, AH[o){Wg) is a Brownian motion with variance (a^ + /3^)A 
per unit time and drift —7, reflected at and H(0). This proves the first 
claim. 

Next, using the definition U^"^^ = $(yV^"^) and the scaling properties of 

it is easy to see that ^?(") = $(>V^"^). Since, by Theorem 3.2, we know 

that yy^, where is continuous and W^(t) has a continuous dis- 

tribution for every t, an application of the continuous mapping theorem, 
along with the continuity property of $ stated in Lemma 4.1, shows that 
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^ This, in particular, implies that C/(") =^(")(M) U* . Since 

[/* has the same distribution as W* , this proves the lemma. □ 

We identify the heavy traffic limit of the workload in the reneging system. 
We start with Proposition 6.2, which states that the number of customers in 
the EDF system with reneging having lead times not greater than the current 
frontier and the work associated with these customers are negligible under 
heavy traffic scaling. Then, in Corollary 6.3, we use the comparison results 
established in Section 5.1 to show that the workloads in the reference and 
reneging systems are equal with high probability and so their heavy traffic 
limits coincide. 

Proposition 6.2. The processes >V(")(0,F(")] and Q(")(0,F(")] con- 
verge in distribution to zero as n— t-oo. 

This result holds for the same reason that state-space collapse occurs 
for priority queues, an idea that can be traced back to [35]. Specifically, 
in our model, due to the nature of the EDF service discipline, the entire 
capacity of the server is always devoted to work that lies to the left or at the 
frontier, as long as the system is nonempty. Thus the process 
is equal to the workload in a single-server GI/G/1 queue that has netput 
process V(")(t)(-oo,F(")(t)] -t, t>0. By showing that < y/Ey* , 

one shows that this (high-priority) queue is in light traffic as n — )■ oo , and 
so its diffusion scaling vanishes in the limit. Since a rigorous proof that 
>V(")[C("),F(")] ^ and Q(") [C^*^), ^ would be very similar to the 
proofs of Proposition 3.6 and Corollary 3.8 in [7], we omit the details. We 

note that VV(")(0,C(")) = Q(")(0,C(")) = by definition. 

Corollary 6.3. Let T > 0. As co, 
(6.8) P[C/(")(t) = VF(")(t),0 < t<nT]^l. 

Proof. Because customers with strictly positive lead times do not re- 
nege, we have W(")('5i"^-)(0, i^^^H-^f ^)) < W(")('Si"^)(0, i^^^n-^i"^)) for k > 
1. Thus, by Lemmas 5.2 and 5.6, to prove (6.8), it suffices to show that as 

n — )• oo, 

P[VV;(")(5W)(o,f(")(5^"))) + vP < 4"\ 1 < A: < A^^'XnT)] ^ 1. 
However, this follows from the fact that, by (2.15), 

max ul"^ = y/n max AN^\t) = o{\/n), 



l<fc<A(")(nr) 0<i<T 



the inequalities > y/ny^, > 0, and Proposition 6.2. □ 



Theorem 3.4 now follows immediately from Lemma 6.1 and Corollary 6.3. 
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6.1.2. Proofs of Proposition 3.5 and Theorem 3.6. We present the proofs 
of the remaining two hmit theorems concerning the measure-valued work- 
load processes. For this, we need two preliminary results. The first, Lemma 
6.4, is that the frontier in the reneging system is strictly positive with high 
probability. The second result. Proposition 6.5, is a recap of a result estab- 
lished in [7]. 

Lemma 6.4. Let T > Q. As n ^ oo, 

(6.9) P[f(") (i) > 0, < t < nT] ^ 1. 

Proof. Let < t < nT. If > 0, then is not smaller than 

the lead time of the currently served customer, so F^'^\t) > 0. If W^"'^{t) = 0, 
then the customer indexed by A^"^^ (t) has already been in service, so 

(6.10) > y/^y^ - n^"^ 



> V^y* — max til"^ 



l<fe<A{")(nT)+l 

However, max^^f,^^(„)^^rps^_^_^u^^^ = o{y/n) by (2.13) (in particular, by the 
fact that S* has continuous sample paths), so (6.10) implies (6.9). □ 

Proposition 6.5 (Proposition 3.4 [7]). Let -co <yo <y* and T > be 
given. As oo, 

sup sup |V(")(t)(y,oo) + i?(y + V^t)-i/(y)| Ao, 

yo<y<y* 0<t<T 

sup sup \A^''\t){y,oo) + XH{y + ^/^t)-XH{y)\-^0. 

yo<y<y' o<t<T 



Proof of Proposition 3.5. Let T > 0. We will show that F^") =^ F* 
in L>R[0,r]. By definition, y* - ^/nt < F(")(t) < y*. Thus, by Proposition 6.5 
and the fact that H{y) = for y > y*, 

(6.11) sup sup |V(")(t)(F(")(t) Vy,oo)-i?(F(")(t) Vy)| Ao. 
0<y<y* 0<t<T 

Putting y = in (6.11) and using Lemma 6.4, we obtain 

(6.12) sup |V(")(t)(F(")(t),oo)-i7(F(")(t))| Ao. 

0<t<T 
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For any t > 0, 

= W(")(t)(0,F(")(t)] + w(")(t)(F(")(t),oo) 
(6-13) _ ^ 

= (t ) (0, (t)] + V(") (t) (t) , oo) , 

where the second Une follows from the fact that none of the customers 
in the EDF system with reneging with lead times at time t greater than 
F^"'\t) has received any service up to time t. This, together with Proposi- 
tion 6.2 and Theorem 3.4, yields V(")(F("), oo) W*. Thus, by (6.12), we 
have H{F^"'^) W* in Dk[0,T]. Applying the continuous function to 
both sides of this relation and using (3.4), we obtain F^"^^ =^ F* in D^[0,T]. 
□ 

Proof of Theorem 3.6. Define a mapping ip:R^ M by the formula 
V'(x)(5) = /Bn[x,oo)(l - <^(^)) dr] for X e M and B G B{R). It is easy to see 
that tp is continuous. Hence, by Proposition 3.5, 

(6.14) (V(F(")), AV(F("))) ^ mF*),XijiF*)) = {W* , Q*). 
Let T > 0. We claim that 

(6.15) sup sup I (i) (y , oo) - V^(f(") (t)) (y , oo) I A 0, 
o<y<y* o<t<T 

(6.16) sup sup |Q(")(t)(?/,oo)-A^(F(")(t))(y,oo)| Ao. 

o<y<y* o<t<T 

Indeed, reasoning as in (6.13), we see that, for < y < y* and < t < T, 
|>vW(t)(y,oo)-i/(FW(t)Vy)| 

< |W(")(t)(F(")(t) Vy,oo) - /f(F(")(t) Vy)| + VV(")(i)(0,F(")(i)] 

= I v(") (t) (f(") (i) V y , oo) - v(f(") (t) ) (y , oo) I + (t) (0, (i)] . 

Therefore, (6.15) follows from (6.11) and Proposition 6.2. A similar argument 
gives (6.16). We have W(")(t)(-oo, 0] = yvW(i)(y*, oo) = Q(")(t)(-oo,0] = 
Q(")(t)(y*,oo) = and, by Lemma 6.4, P[V'(F(") (t))(-oo, 0] = 0,0 < t < 
T] — 7- 1 as n — )■ oo. Also, ^{x){y* ,oo) = for every x E M. Thus, (6.14)~(6.16) 
imply that ()?("), Q^")) ^ (W*, Q*) in Dm[0,T]. □ 

6.2. The heavy traffic limit of the reneged work process. In this section, 
we identify the limit of the sequence {R^^ , n G N}, thereby proving Theorem 
3.8. To do this, it is convenient to show that many of the processes under 
consideration can be put on a common probability space so that certain 
weak limits established earlier can be replaced by almost sure limits. 
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Lemma 6.6. The processes yV^^\ n E N, yV*s, U* and W* 

can be defined on a common probability space {Q,J-,¥) such that P almost 
surely, as ?i — )• oo, 



(6.17) wP^W*s, 

(6.18) W^^'^^W^, 

(6.19) wP(.)(^-^,o]^w*s{-){-^,0] = mi-)-Hm^, 

(6.20) f/^") ^ U* 
and 

(6.21) |^(")^VF* = ;7*, 



where W^"-^ = W^iR), = >V^(M), C/W =^?(")(M) and U* =U*{R). Fur- 
thermore, Wg is a Brownian motion with variance (q^ + /3^)A per unit time 
and drift —7, reflected at 0, while U* is a doubly reflected Brownian motion 
on [0,H{0)], also with variance (a^ + /3^)A per unit time and drift —7. In 
particular, 

(6.22) U* = AHio){W*s) = W*s-KX + K*_, 

where K'^. are the unique RCLL nondecreasing functions satisfying ^^'^(O) = 
and 

(6.23) / I{U'is)<Hio)}dK*+{s) = oJ I{U'is)>o}dK*_{s) = 0. 

J[0,oo) J[0,oo) 

The almost sure limits in (6.17)~(6.21) hold uniformly on compact intervals. 
Proof. Recall from Theorem 3.2 that W^. Using the Sko- 

In) (n) 

rokhod representation theorem, we construct the model primitives , Vj 
and L^"^ for j G N and n € N on a common probability space (ri, P) such 

' — 'in) 

that the sequence of processes , n G N, and the limiting process 
are defined on this space and (6.17) holds. Here and below the almost sure 
convergences are in the Ji topology on D_m [0, 00) or -Dm[0, 00), and since the 
limits are continuous in every case, this is equivalent to uniform convergence 
on compact intervals. Since the mapping / : [0, 00) 1— )■ -De[0, 00) given by 
= m(')(^) is continuous, we have (6.18). Under P the measure- valued 
process constructed on has the same distribution as the process 
appearing in Theorem 3.2, and thus takes values in the set of measure- 
valued process of the form jBn[F°it) oa)^'^ ~ ^iy))dy for some RCLL pro- 
cess F°(i). However, T^*(t) = /iRnVs(t),oo) " Giy))du = HiF^it)); hence 
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F^{t) = F^{t) is given by (3.2). In other words, with defined by (3.2), 
the first equation in (3.3) holds. Due to Proposition 3.1, the above argument 
also shows that under P, is a Brownian motion with variance (a^ + /3^)A 
per unit time and drift —7. In addition, since for each t, the measure Wg{t) 
is nonatomic, we have (6.19). 

Now, following (4.1) and (6.3), we set ZY^") = $(wj*^) and U* = $(W^). 
Also, as defined in (6.2), let U^"'^ be the scaled version of^/("), and let f/(") 
and U* be as defined in the statement of the lemma. Then l/{^"'\U^"'\n G N, 
U* and U* are also defined on {Q,J^,¥) and (6.20) follows from Lemma 4.1. 
This implies (6.21). Since U* = $(>V^)(R) = Aj^(o)(W^), the characteriza- 
tion of U* as a doubly refiected Brownian motion that satisfies relations 
(6.22) and (6.23) is a consequence of the statements following (6.4), in par- 
ticular, (6.5) and (6.6). 

Since the model primitives u^"^ , and L^""^ for j G N and n G N are all 
defined on {fl,T,F), so are the workload process W^^'^ and its scaled version 
I^("). Corollary 6.3 implies that and have the same limit, and 

hence (6.21), the almost sure counterpart to Theorem 3.4, holds. □ 

The assertion of Theorem 3.8 is that 

(6.24) R^w^^K^ 

where K"^ is the local time for U* at H{0) from (6.22). For T < 00, define 

(6.25) Zn{T) ^ 0^\t) = R\;^\t),0 <t<T}. 

Prom the workload evolution equations (4.69) and (4.70), it follows that if 
t/(")(t) = W?W(t) for t G [0,r], then R'l^\t) = ^[^^(t) for t G [0,r]. Hence, 
by Corollary 6.3, we know that for every T < 00, ¥{Zn{T)) — )• 1 as ?i — )■ 00, 

which shows that the limits in distribution of R^ and R\y , n G N, must 
coincide (if they exist). Purther, since 

^1") = R^^^ by Corollary 4.9, these 

must be equal to the limit in distribution of , n G N. Hence, to complete 
the proof of Theorem 3.8, it suffices to show that 

(6.26) K^^^^Kl. 

For n G N and k>l, recall the definitions of t^^\ and given in (4.11) 

and (4.12), respectively, and define t^'L^ = ^t^^i and ct^"^ = ^c^"^- Applying 
the heavy traffic scaling to (4.13), it is easy to see that for t>0, 

(6.27) 



Ki^\t) = Y,m-\a'i'^-)V max wj^^ (.)(-oo, 0] 
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Keeping in mind the limits in (6.17) and (6.19), we introduce the related 
process 

(6.28) yW(t)^j;[T^|(ai"V max (W*s{s) - Hm"- -W*s{dP] 

for t > 0, and denote the difference by 

(6.29) e(") (i) = y (t) - i?!"^ (i) Vt > 0. 

Then is nondecreasing and continuous, and e*-") is an RCLL process. 

In the next two lemmas, we show that y^"-) increases only when U* is 
at H(0) and that the difference e^^'^ between y("') and is negligible in 
heavy traffic. The main reason for introducing the sequence yW, neN, IS 
that it facilitates the proof of the former property. 

Lemma 6.7. For every n G N, y('") and are constant on each in- 

terval [t^i^}-^ , S"^""* ) , k>l. Moreover, 

(6.30) / V(,)<H(o)}dy(")(i) = 0. 

J[0,T] 

Proof. Fix n G N. The first statement follows immediately from (6.27), 
(6.28), and the fact that the intervals [T[."_^-|^,a^"^) and [a^"\r^"^), A; > 1, 
form a disjoint covering of [0,oo). Now, fix fc > 1 and let be the set of 
points t G [S'^"'')^^!"^) such that 

(6.31) W*s(a'-^^)< max (W^(s) - i/(0))+ = M^|(i) - if(0). 

Since Wg is continuous, J^"^ is closed, and so its complement in ''"^"^ ) 
is the union of a countable number of open intervals, with possibly one 
half-open interval of the form [a^^^a) for some a > 5"^"^ From the explicit 
formula for y(") given in (6.28), it is easy to deduce that y^") is also constant 
on each such interval. Thus, to establish (6.30), it only remains to show that 
for each k> 1, 



(6.32) / huHt)<Hm}dY(^Ht) = 0. 



Fix t G J^"^ and note that by the equality in (6.22) and the definition 
(6.4) of Ah{o), we have U*{t) = W^{t) - K*{t), where 

(6.33) K*{t)= sup \{Ws{s) - H{0))~^ A inf Ws{u) 

se[0,t] L u£ls,t] 
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Also, note that 




Since K* (t) is equal to the maximum of the quantities on the left-hand side 
of the last two displays, we conclude that 

K*it) < W*siaP) V {W*sit) - H{0)) = W*s{t) - i/(0), 

where the equality follows from the inequality in (6.31). This, when combined 
with the fact that U*{t) G [0,H{0)], shows that U*{t) = W^{t) - K*{t) = 

H{0) for all t G which proves (6.32). □ 

We recall some standard definitions that will be used in the next lemma. 
Given / G I^[0, oo) and < ti < ^2 < oo, the oscillation of f over [ti,t2] is 



Proof. Fix T > and let 77 > be arbitrarily small. By the Kolmogorov- 
Centsov theorem (see, e.g.. Theorem 2.8, page 53 of [17]), we can construct a 
positive, increasing deterministic function 6{-) satisfying liuisio 0(5) = and 
majorizing the modulus of continuity ww* (s [Oj ^]) of the reflected Brownian 
motion W* over [0,r] on a set J7 with P(0) > 1 - r?. 

For each subsequence in N, there is a sub-subsequence S along which the 
limits (6.17)-(6.21) hold P-almost surely. We choose so that these limits 
hold along S for all w G fi. 

In what follows, for n G 5, we denote Zn{T) simply by Zn, and evaluate 
all processes below at a fixed w G n il. Choose A < ?/=k/3, and let no G 5 



Osc(/; [tiM) = sup{|/(t) - f{s)\:ti <s<t< ts}. 



and the modulus of continuity of f over [0,T] is 



Wf{6; [0,r]) ^ sup{|/(t) - f{s)\ :0<s<t<T,\t-s\<6}. 




p 
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be such that for all n £ S, n > no, 

(6.34) sup \wP{t){-oo,0] - {Wsit) - H{0))^\ < A, 

t£[0,T] 

sup \wi!'\t-)-wm<^, 

te[o,T] 

(6.35) 

sup \W^''\t-)-W*{t)\<A, 
te[o,T] 

(6.36) sup |[/(")(i-) - U*{t)\ < A. 

t6[0,T] 

From the definitions (6.27) and (6.28), respectively, of and y^") it is 

(n) 

clear that, for every /c G N such that < T, 

sup |y(")(i) - i^^"H4"^-) - 0+\t) - K^+\dP-))\ < 2A. 

Define 

J„ ^ {^^ G N:^f )(f(")) - ^i")(5l")-) > 0,?(") < T}, 
J„ = {A: G N:y(")(?l")) - y(")(at")-) > 0,?!"^ < T}, 

and let c^") be the cardinality of U J("). Since and y(") are both 

constant on intervals of the form [t^'I^ , ^[.'^^ ) , k>l (see Lemma 6.7), we 
have 

(6.37) e('^)(r) = sup |y('^)(s) - ^i"^(s)| < 2c(")A. 

se[o,T] 

We now claim that 

(6.38) ke[JnUJn] Osc(VF*,[ai"\fi")])>|l. 

We defer the proof of the claim and instead first show that the lemma fol- 
lows from this claim. Let 9~^{-) denote the inverse of 9 and define M = 
T/9~^{y^:/3) < oo. From the claim, we conclude that if A: G [Jn U J„] then 

"^fc"^ - ^1"^ ^ 0~^{y*/3) > 0, which in turn implies that c^") < M. Substitut- 
ing this into (6.37), we conclude that for every A > 0, there exists no(A) G 5 
such that for all n G 5, n > no (A), 

P(e(")(r) > 2MA) < F{Z^ U 0") < P(Z^) + r/. 

Taking limits as n — )■ oo through S and using the fact that ¥(Zn) — )• 1, 
we conclude that e*-"'(T) — )■ 0. We have shown that for each subsequence 
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in N, there is a sub-subsequence along which e*-"'(T) — > 0. It follow that 

£i'^){T) — ^ 0, where the limit is taken over all n G N, and this proves the 
lemma. 

We now turn to the proof of the claim (6.38). Note first that by the defi- 
nition of -ff(O) and y^,, we have -ff(O) > y*. If A; G J„, then Lemma 6.7 shows 
that U*{t) = H{0) for some t G [5^"^?^"^). By the equality = 

proved in Lemma 4.7 and (6.36), this implies that the oscillation of U* on 

[^fc'^^i"'*) ^^^^ than H{0) — A > y*/3. Since W* = U*, the conclusion 
in (6.38) holds. 

Finally, suppose k G Jn- Since ^ — ^ — -^vi/\ have 

that is, the deadline of a customer in the reneging system expires during the 
unsealed time interval Since 

(6.39) Ty(")(4")-) = 

[because U'^^\a'^^^ —) = and, by Lemma 5.2, < C/^")], this customer 

must arrive during the interval [c^^^ r^""* ) . Since his initial lead time is 

greater than or equal to y^y* , there is a time n^o £ [^i"^ ) '^i"^ ) when this 
customer has lead time exactly -y/ny*. After time nig, this customer cannot 
be preempted by new arrivals, all of which have initial lead times greater 
than or equal to ^/ny^. At time nto, the work that must be completed before 
this customer is served to completion is at most W^"^(nio)(0, -^/Hy*]. Since 
this customer becomes late, we must have VF(nto) > VV''"^(nto)(0, -^/ny^,] > 
i/ny*, or equivalently, VF^")(to) > VV*-"'^(to)(0,y*] > y*- By right continuity, 
W^^\{tQ + y)—) > y* for some > so small that to + v < t^j^^ . From the 

second inequality in (6.35) and the fact that W^"'^ (a^j^^ —) = [the scaled 
version of (6.39)], we conclude that 

w*{to + u)-w*{^^^^)>^, 

and this gives us the conclusion in (6.38). □ 

Proof of Theorem 3.8. Fix T < oo. Let (^(") =U* - C/("), and let 
=sup,g[o,T]|t^*(s) -f^^^Hs)!- According to (4.6) and (4.15), 

We scale this equation to obtain 

(6.40) U* = W^"'^ + (^('^^ - ^i") + Ki"^ = t?^"^ + 6^^'^ + - Y^^'^ + K^^\ 
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where [cf. (4.14)] 



km 



K^"' is defined by (6.27) and e^") is defined by (6.29). According to (4.16), 

fo ^{c7(")(i)>o} 



fo'^(u(")(t)>o\ dK'f!\t) = 0, which imphes 



(6.41) [\u^,t)>t-^}'^-'^') = ^- 

Since + J^") + e^") W*g due to (6.18), (6.20) and Lemma 6.8, and, 
by (6.22), U* is obtained by applying the Skorokhod map on [0,if(0)] to 
VF^, the convergence (6.26) is an immediate consequence of (6.40), (6.41), 
Lemmas 6.7, 6.8 and the invariance principle for reflected Brownian motions. 
However, since we are in a particularly simple setting here, we will provide 
a direct proof without invoking the general invariance principle. 

We choose no so that (^("o) < H{Qi) /3 and recursively define stopping times 
po = 0, and for A: > 1, 



Vk = va.va.< t > pk-i 



3 



3 



Then = po < i^i < pi < ^2 < ■ ■ ■ and linifc^^oo Pk = hmfc^oo i^k = oo. For n > 
no, -ftri"^ is constant on each of the intervals [i^k^Pk]- Moreover, Lemma 6.7 
implies that for each k, y('^) is constant on each of the intervals [pk-ij'^k]- 
For t£ [i^k,Pk], we have from (6.40), (6.18), (6.20) and Lemma 6.8 that 

y(")(t)_i>W(j.^) = w^"^)(^t)-U*{t)+6^''\t) + e^''\t) 

- W^^'H^k) + U*{uu) - S^'^Xyk) - e^"H''k) 

A Wm - U*{t) - {W*s{vu) - U*{vu)). 
It follows that, uniformly for t G [0,T], y(")(t) converges in probability to 

(6.42) Y.^W*s{{t V Uk) A Pk) - U*{{t V i/fc) A pk) - {W*sii^k) - U*iuk))]. 

km 

However, (6.23) implies that for each k, is constant on [ui^jPf^], and K"^ 
is constant on [pk-ij^^k]- Therefore, (6.22) implies that for t £ [/^/t,/Ofc], 

Klit) - Klivk) = W*s{t) - U*{t) - {W*s{uk) - U*{uk)). 

This implies that the expression in (6.42) is K'^(t). But and K^^'^ have 
the same limit in probability because of Lemma 6.8, and we conclude that 

(6.43) max \Ki\t) - Kl(t)| A 0. 
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Convergence in probability implies weak convergence, and we have (6.26). 
□ 



7. Performance evaluation and simulations. We use the heavy traffic 
approximations of this paper to evaluate the performance of the system with 
reneging and compare this to the system in which all customers are served 
to completion. The predictions of the theory, derived in Section 7.1 and 
compared to simulations in Section 7.2, are predicated on the assumption 
that one can interchange the limit as n — )■ oo and the limit as time goes 
to infinity of the fraction of reneged work. A formal proof would require a 
coupling argument such as that found in [37]. The simulation results attest 
to the accuracy of the approximations derived in Section 7.1 and also show 
the great difference in performance between the reneging and nonreneging 
systems. 



7.1. Derivation of theory predictions. We derive formulas (1.1)-(1.7). We 
begin with one of the main results of this paper. Theorem 3.4, which states 
that the limiting scaled workload in the reneging system is a reflected Brow- 
nian motion in [0, //(O)] with drift. More specifically, 

(7.1) W*{t) = W^{t)- Kl{t) + K*_{t), 

where Wg{t) is a refiected Brownian motion on [0,oo) with variance cr^ = 
A(a^ + /3^) per unit time and drift —7, -ftT* is the nondecreasing process 
starting at i^'* (0) = that grows only when W* = 0, and is the nonde- 
creasing process starting at K'^{0) = that grows only when W* = H{0). 
We further saw in Theorem 3.8 that K1^{t) is the limit of the scaled work- 
load that reneges prior to time t in the diff'usion scaling, that is, ^/nK1^{t) is 
approximately the (unsealed) workload that reneges in the nth system prior 
to time nt. 



Lemma 7.1 ([15], Proposition 5, page 90). We have 

7 



(7.2) lim -Kl{t) = lim -EKl{t) = { 

t—^oo t t—^co t 



C72 



2H{0) 



if 1^0, 
z/7 = 0. 



Proof. The first equality in (7.2) is a consequence of the fact that 
W* has a stationary distribution [see (7.5) below]. For the proof of the 
second equality, recall that Wg has the decomposition (2.16). Let / be a 
function. Applying Ito's formula to f{W*{t)) and taking expectations, we 
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obtain 

fXO)E[rsit) + K*_it)] - nHiO))EKl{t) 



(7.3) =E 

+ EfiW*{t))-fiO). 



lf'iW*is))-^a'f"{W*is)) 



ds 



Taking f{x) = x,we obtain E[rg{t) + K*_{t)] -EKl{t) = -ft + EW*{t) - /(O). 
If 7 7^ 0, we may take f{x) = ^e^'^^/'^'^ in (7.3), which leads to the equation 
E[I*{t) + Kl{t)] - e2^-f^(0)/'^'EK;(t) = |i(Ee2^^*W/'^' - 1). Solving these 
equations for EK'^{t), we obtain the second equality in (7.2) for 7 / 0. To 
obtain this equality for 7 = 0, we take f{x) = x^ . □ 

According to (2.12), the work that arrives to the nth system by time nt is 
= ^iV(")(t) + nt. But, iV(") is approximately N* , and hence 

JnN'<^\t)^nt , JnN*(t)^nt ( 7 
lim w lim — ' 



t-^co nt t~^oo nt \ , 

Therefore, if 7 7^ 0, the long-run fraction of reneged work is approximately 

t^oo y{«)(^(")(nt)) ^t^oot t^oo\ nt J 



fn 



(l-7/V^)(e27^^(0)/a2 

Finally, (2.4) implies that the expected lead time in the nth system is 
Ehf^ = /q°°(1 - G{y/y/n)) dy = V^i?(0). Using this formula and (2.10), we 
conclude that the fraction of work that reneges in the nth system when 7 7^ 
is approximately 



(7.4) 
where 



1 - 1 - p(") 



(„)^g2(i-p{-))E4"V-2 _ 1) p(")(e«^ - 1) ' 



P' 



We have suppressed the dependence of Q and D on n, which will remain 
fixed. If 7 ■ 
and (1.2). 



2 

fixed. If 7 = 0, then in place of (7.4) we have We have established (1.1 
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Remark 7.2. Corollary 3.7 also implies that the limiting scaled queue 
length process is ATF*, which is a doubly reflected Brownian motion in 
[0, Aff(O)] with drift — 7A and variance per unit time A^o"^. This incorrectly 
suggests that X^/nK^{t) is approximately the number of customers who 
renege in the nth system prior to nt. The simulations indicate that this 
naive interpretation of Corollary 3.7 applied to the queue length process is 
incorrect, as does the following heuristic. 

According to [15], Proposition 5, page 90, if 7 7^ 0, the stationary density 
for W* is 

{2^g-27x/a2 
— — ^ , . , if 0<a;<iJ(0), 
0, otherwise, 

whereas the stationary density is uniform on [0, H{0)] if 7 = 0. Therefore, 
for 7 7^ and t large, the density of W^''^\nt) w ^/nW*{t) is approximately 

( 6e~^^ — 
otherwise. 



n 



We have suppressed the dependence of (/? on n. 

Suppose now that the lead times of arriving customers are not random. 
Then in the nth system, all lead times are equal to y/nH{{)) = D. In this 
case, the EDF policy serves customers in order of arrival (FIFO). Suppose 
the workload in queue is W at the time of arrival of a customer whose 
service requirement is V. Recall that the expected service time is 
and because n is fixed, we suppress it and write = The arriving 
customer will be served to completion if and only \i W + V < D. Suppose 
further that the arrival process ^^"^ is Poisson, so that according to the 
PASTA property ("Poisson arrivals see time averages"; see [1], Theorem 
6.7, page 218), an arriving customer will encounter a workload W having 
approximately the distribution (p. The probability the arriving customer 
eventually reneges is thus 



"{W >D-V} = ¥.[¥{W >D-V\V}]=¥. 



D 



_ f{w) dw 
(D-V)+ 



Because D is of order ^/n and V is of order 1, we have {D — V)~^ = D — V 
with high probability. Using this approximation, we complete the calculation 
for the case 7 / to obtain 

(7.6) P{Customer reneges} ^ (Ee^^ - 1). 

g6»D _ I 
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If the customer reneges, then work V + W — D > is lost. The expected 
lost work is 



E.[V + W — D I Customer reneges] w E 



I$-V)+ (y + w- D)ip{w) dw 



P{ Customer reneges} 
Again using the approximation (D — V)^ ^ D — V, we obtain 

E[F + — D I Customer reneges] 



1 KV 

~ Ee^^ - 1 

1 EV 



9 eEV+l/2e^E[V^] + 0{n-'^/^) 

_ E[V^] 
2EV ■ 

The last expression is, perhaps not surprisingly, the formula for the average 
residual lifetime of a renewal cycle (see [32], Example 3.6(b), pages 80 and 
81). Consequently, when lead times are constant and the arrival process is 
Poisson, we should expect the total number of customers reneging in [0, t] 
times the expected amount of work lost per reneging customer to approxi- 
mately equal the total amount of work lost by reneging in [0,t]. If we divide 
both by the total number of customer arrivals in [0,t] and take limits as 
t — )• cxD , we find 

Fraction of lost customers in reneging system 

Fraction of lost work in reneging system „ 

(7.7) « = *EV 

E[V + W — Z) I Customer reneges] 

2(Ey)2 

~ — - — 7T— X (Fraction of lost work in reneging system) . 
E[\/ ^] 

This is (1.7) with EV = ^ and E[V'^] =f3^ + jj. 

If y is exponentially distributed, hence E[y^] = 2{EV)'^, then (7.7) implies 
that the fraction of customers who renege will be approximately the fraction 
of work that reneges. See Figure 1 for simulations that confirm this assertion. 
On the other hand, if V is nonrandom, hence equal to its mean l//i, then 
(7.7) predicts that the fraction of customers who renege will be twice the 
fraction of work that reneges. See Figure 2 for simulations that confirm this 
assertion. Both these conclusions hold irrespective of the value of A. 

The last conclusion is inconsistent with a naive interpretation of Corollary 
3.7, according to which work reneges at a rate 1/A times the rate of customer 
reneging. Since work arrives at a rate EV ~ 1/A times the rate of customer 
arrivals, this naive interpretation of Corollary 3.7 would say that the fraction 
of work reneging would approximately agree with the fraction of customers 
reneging regardless of the distribution of V. 
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Fig. 1. M/M/l queue. 
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Fig. 2. M/D/1 queue. 
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We next turn our attention to the performance of the standard (nonreneg- 
ing) system. Recah from (2.15) that the scaled workload process when all 
customers are served to completion converges to Wg, a reflected Brownian 
motion with drift —7 (we now assume 7 > in order to have a station- 
ary distribution) and variance o"^. In particular, Wg"'\nt) ~ ^/nWg{t). The 
stationary density for Wg is 

*f .A j ^e-^W'^', if x>0, 
1 0, otherwise, 
and so for large t, the density of W^^\nt) is approximately 

^ ' 10, otherwise. 



Consequently, the long-run fraction of time W^^^ spends above level D is 
g-SD rpj-^g workload level at which the limiting frontier reaches is H{0), 
and hence it is approximately the case that the nth system sees lateness if 
and only if T^(") exceeds D = y/nH{0). In other words, the theory predicts 
that 

Fraction of late customers in standard system 
(7.8) = Fraction of late work in standard system 

We are using here the result for GI/G/1 queues that 



lim lim — / I dt 

1 

= lim lim — / I,rj;(„),^, 



1 



T 







a result that grows out of the work of Kingman [19, 20] (see [10] for a general 
result that specializes to the case under consideration). 

It is important to compare the fraction of work that reneges in the reneg- 
ing system, given by (7.4), with the fraction of work that is late in the 
standard (nonreneging) system. The ratio of these quantities of lost/late 
work is 

Lost work in reneging system e^^ { ^ ~ P^"''^ 



Late work in standard system gSD — \\ pi"^) 
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The parameter is 0(1/VH), OD is 0(1), and 1 - is 0(1/^). Thus 
the ratio in (7.9) is 0{\/^/n). 

Remark 7.3. If lead times are a nonrandom constant EDF reduces 
to fir st-in- first-out, and the fraction of lost customers in an M/G/1 queue 
with < p < 1 is (1 - p)^{W > D}/{1 - pP{W > D}), where W is the 
steady-state workload in the corresponding nonreneging M/G/1 queue (see 
[3]). In the heavy traffic limit of our model, ¥{W > D} = e~^^ [see the 
derivation of (7.8)]. Recalling that 1 — p = 0{1/ ^/n) in (1.1), we observe 
that this is consistent with (1.1). 

7.2. Simulation results. We conducted a simulation study to assess the 
accuracy of these approximations and to compare the performance of the sys- 
tems with and without reneging. Two systems were considered, an M/M/1 
system presented in Figure 1 and an M/D/1 system presented in Figure 2. 
In both cases, A = 0.5 and ^ = 1.96, and so the traffic intensity is p = 0.98. 
These parameter values result in 6 = 0.010202 for the M/M/1 case and 
= 0.02 for the M/D/1 case. The initial deadline distribution is uniform on 
[5, B] with the mean deadline D = varying from B = 5 (constant dead- 
lines) to i? = 200. The data points are the simulation results averaged over 
one billion customer arrivals per case. The curves that are superimposed on 
the data are the theoretical values, e~^^ for the case in which customers are 
served to completion (the standard system), and equations (1.1) and (1.7) 
for the fraction of work lost and the fraction of lost customers for the reneg- 
ing system. Equation (1.7) is derived in Remark 7.2 under the assumption 
of constant deadlines. Nevertheless, we apply it for the variable deadline 
case in the simulation study. The fraction of late work or late customers for 
the system in which customers are served to completion is also presented to 
compare its performance with that of the reneging system. 

The M/M/1 results are presented in Figure 1 with the fraction of cus- 
tomers missing their deadlines, the fraction of customers reneging, and the 
fraction of work reneging plotted on a log scale on the y-axis against the 
mean deadline on the x-axis. There is nearly perfect agreement between 
the theoretical approximation and the simulation. In fact, one cannot see 
the plot of "Fraction of Customers Late (No Reneging)" because it coin- 
cides with the "Theory" plot at the top of the figure. Similarly, one can see 
only parts of the plots of "Fraction of Customers Reneging" and "Fraction 
of Work Reneging" because they coincide with the "Theory" plot in the 
middle of the figure. One can see the linear form for the case of service to 
completion. Furthermore, the simulation confirms the prediction of (1.1)^ 
(1.4) that for sufficiently large values of D, the performance of the reneging 
system is parallel on a log scale to that of the standard system with the two 
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curves separated by approximately 0.02. This corresponds to a reduction in 
work that misses its deadline by a factor of 40 to 50. 

Figure 2 presents the results for the M/D/1 system. The results are qual- 
itatively identical to those of Figure 1, except the fits of the theoretical 
curves are not as exact as the fits for the M/M/1 system; it appears that 
now the value = 0.02 is slightly too small and hence the theory slightly 
overestimates the fraction of work that misses its deadline, especially when 
the mean deadline is large. Also, the lost or late work and the customer loss 
or lateness fractions are significantly smaller than for the M/M/1 system 
owing to the reduction in variability of the customer service time distribu- 
tion. The reduction in missed deadlines between the two systems for large 
values of D is again a factor of 40 to 50. In both figures, it is clear that 
there are one to two orders of magnitude of improvement in the overall per- 
formance of the system resulting from stopping service on customers when 
their deadlines expire. 

APPENDIX: OPTIMALITY OF EDF 

Proof of Theorem 5.1. Let vr be a service policy and let to be the 
first time vr deviates from the EDF policy, either because it idles when there 
is work present, or it serves a customer other than the customer present with 
the smallest lead time. Let j be the index of the customer with the smallest 
lead time at time Iq. 

We consider first the case that vr idles at time to- Ii^ this case, we define 
p(vr) to be the policy that emulates vr except as noted below. From time 
to, whenever vr idles, /^(vr) serves customer j, at least until time ti, when 
customer j leaves the /9(vr) system because either /9(vr) serves customer j to 
completion or else the deadline of customer j elapses. From time ti, p(vr) 
idles if vr serves customer j. We will show that for t > 0, 

(A.l) i?^(,)(t)<i2^(t). 

Let Vk{t) [resp., f^(t)] be the residual service time of the kih. customer at 
time t under vr [resp., /o(vr)]. In particular, if is the deadline of the kth 
customer, then Vk{dk — ) [resp., — )] is the work corresponding to this 

customer that is deleted by vr [resp., /3(vr)] due to lateness, and 

(A.2) i2pw(t)= R^{t)= J2 Mdk-)- 

k: d^.<t k: d]^<t 

By the definition of p(vr), for t > and k^ j, we have 
(A.3) v'^{t) = vk{t), 

whereas 



(A.4) 



v''it)<v,it). 
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Summing (A. 3) over k^j, invoking (A. 4) and (A. 2), we obtain (A.l). 

We next consider the case that at time to, vr serves customer j- In 
this case, we define p{Tr) to be the policy that emulates vr except as noted 
below. From time to, whenever vr serves customer i, p{Tr) serves customer 
j, at least until time ti, when p{'jr) serves customer j to completion or the 
deadline of customer j elapses. Prom time ti, p{7r) serves customer z if vr 
serves customer j, provided customer i is present in the system under p{tt). If 
vr serves customer j and customer i is not present under p{7r), then p{7r) idles. 
We again have (A. 2) and (A. 4), whereas (A. 3) now holds only for k ^ 
If the ith customer is served to completion under p(vr), then vf{di — ) = 0, 
and (A. 3) for k ^ and (A. 4) imply that (A.l) holds for all t. It remains 

to consider the case that the ith customer becomes late under p('7r). In this 
case (A.3) for k ^ {i,j} and (A.4) imply that (A.l) holds for t € [0,di). Let 
wi denote the work done by p{tt) on the jth customer when vr works on the 
ith customer in the interval [to,ti). Let W2 be the work done by p{n) on 
customer i in the time interval [ti,oo) while vr works on customer j in this 
time interval. Finally, let be the work done by tt on customer j in the time 
interval [ti,oo) while p{tt) is idle. Then Vj{dj—) + wi = Vj{dj—) +1^2 + ^3 
and v1{di—) + W2 = Vi{di—) + wi, which implies 

(A.5) v''^{dj-) + v1{di-) = Vj{dj-)+Vi{di-) + W3. 

We argue by contradiction that cannot be positive. If were positive, 
then at some time t>ti, tt serves customer j and customer i is not in the 
p{tt) system. This implies that dj > t, and since by assumption, di > dj, the 
absence of customer i in the p{tt) system means that this system has served 
customer i to completion. We conclude that v^{di—) = 0. On the other hand, 
customer j is also not in the p{Tr) system at time t>ti, and so Vj{dj—) = 
as well. The left-hand side of (A.5) is zero, and hence W3 must be zero. We 
conclude that 

(A.6) vP{dj-)+vP{d^-) = Vj{dj-) + Vi{d,-). 

Since dj < di, if t > di, then (A.3) for k ^ {i,j} and (A.6) imply (A.l). 

Starting from the service policy tt, we have obtained a service policy p{tt) 
that either is work conserving until the departure of customer j or else 
gives customer j priority over customer i until the departure of customer j. 
However, immediately after time to, the policy vr may serve some customer 
k ^ {i, j}, and hence p{tt) also serves k at this time, although customer j 
is more urgent. Therefore, we apply n iterations of the mapping p, where n 
is the number of customers in the tt system at time to, and thereby obtain 
a policy that is work-conserving and serves in EDF order at least until the 
first time after to that there is a departure or an arrival. We have Rpn(^^'j{t) < 
R^{t) for ah t > 0. 
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By assumption, for each t the number of system arrivals A{t) by time t 
is finite. Hence the maximum number of customers in the system over the 
interval [0,t] is bounded by A{t), and the number of arrivals and departures 
up to time t is bounded by 2A{t), irrespective of the service policy. Thus, 
if we start with any policy vr, the number of iterations of the mapping p 
required to obtain a policy that is work conserving and serves in EDF order 
up to time t is finite. Under this policy the amount of work removed by 
lateness up to time t is the same as for the EDF system in the theorem, and 
hence (5.1) holds. □ 

Remark A.l. In the above proof we have implicitly assumed that vr 
[and thus /o(7r)] never serves more than one customer at the same time. This 
assumption simplifies the exposition of the argument, and the generality of 
Theorem 5.1 is sufficient for this paper. However, the proof can be general- 
ized to policies permitting simultaneous service of customers (e.g., processor 
sharing). In this case, in the construction of p{Tr) we must additionally take 
the rates at which customers receive service into account. For example, the 
difference in the rates with which the jth customer receives service under 
p(7r) and tt in the time interval [to,ti) must be equal to the rate of service 
of the ith. customer under tt in this time interval, the rates of service of all 
other customers in this time interval under vr and p('/r) must be the same, 
etc. 

REFERENCES 

[1] ASMUSSEN, S. (2003). Applied Probability and Queues, 2nd ed. Applications of Math- 
ematics (New York) 51. Springer, New York. MR1978607 

[2] BiLLlNGSLEY, P. (1999). Convergence of Probability Measures, 2nd ed. Wiley, New 
York. MR1700749 

[3] Boots, N. and Tijms, H. (1999). A multiserver queueing system with impatient 

customers. Management Science 45 444-448. 
[4] Chen, H. and Mandelbaum, A. (1991). Stochastic discrete flow networks: Diffusion 

approximations and bottlenecks. Ann. Probab. 19 1463-1519. MR1127712 
[5] Decreusefond, L. and Moyal, P. (2008). Fluid limit of a heavily loaded EDF 

queue with impatient customers. Markov Process. Related Fields 14 131-158. 

MR2433299 

[6] Down, D., Gromoll, H. C. and Puha, A. (2009). Fluid limits for shortest remaining 
processing time queues. Math. Operations Research 4 880-911. 

[7] DOYTCHINOV, B., Lehoczky, J. and Shreve, S. (2001). Real-time queues in heavy 
traffic with earliest-deadline-first queue discipline. Ann. Appl. Probab. 11 332- 
378. MR1843049 

[8] DuPUlS, P. and Ramanan, K. (1999). Convex duality and the Skorokhod problem. 

Probab. Theory Related Fields 115 197-236. 
[9] Ethier, S. N. and Kurtz, T. G. (1986). Markov Processes: Characterization and 

Convergence. Wiley, New York. MR838085 



64 



KRUK, LEHOCZKY, RAMANAN AND SHREVE 



[10] Gamarnik, D. and Zeevi, A. (2006). Validity of heavy traffic steady-state ap- 
proximation in generalized Jackson networks. Ann. Appl. Probab. 16 56-90. 
MR2209336 

[11] Gromoll, H. C. (2004). Diffusion approximation for a processor sliaring queue in 

fieavy traffic. Ann. Appl. Probab. 14 555-611. MR2052895 
[12] Gromoll, H. C. and Kruk, L. (2007). Heavy traffic limit for a processor sharing 

queue with soft deadlines. Ann. Appl. Probab. 17 1049-1101. MR2326240 
[13] Gromoll, H. C., Kruk, L. and Puha, A. Diffusion limits for shortest remaining 

processing time queues. Preprint, Dept. Mathematics, Univ. Virginia. Available 

at http : //arxiv . org/abs/1005 . 1035 . 
[14] Harrison, J. M. and Reiman, M. Reflected Brownian motion in an orthant. Ann. 

Probab. 9 302-308. 

[15] Harrison, J. M. (1985). Brownian Motion and Stochastic Flow Systems. Wiley, New 
York. MR798279 

[16] Iglehart, D. L. and Whitt, W. (1970). Multiple channel queues in heavy traffic. 

I. Adv. m Appl. Probab. 2 150-177. MR0266331 
[17] Karatzas, I. and Shreve, S. E. (1988). Brownian Motion and Stochastic Calculus. 

Graduate Texts in Mathematics 113. Springer, New York. MR917065 
[18] Kaspi, H. and Ramanan, K. (2011). Law of large numbers limits for many-server 

queues. Ann. Appl. Probab. 21 33-114. 
[19] Kingman, J. F. C. (1961). The single server queue in heavy traffic. Proc. Cambridge 

Philos. Soc. 57 902-904. MR0131298 
[20] Kingman, J. F. C. (1962). On queues in heavy traffic. J. Roy. Statist. Soc. Ser. B 

24 383-392. MR0148146 
[21] Kruk, L. (2007). Diffusion approximation for a G/G/1 EDF queue with unbounded 

lead times. Ann. Univ. Mariae Curie- Sklodowska Math. A 61 51-90. 
[22] Kruk, L., Lehoczky, J. P. and Shreve, S. (2003). Second order approximation 

for the customer time in queue distribution under the FIFO service discipline. 

Ann. Umv. Manae Cune-Sklodowska Sect. AI Inform. 1 37-48. MR2254131 
[23] Kruk, L., Lehoczky, J. P., Shreve, S. and Yeung, S.-N. (2004). Earliest-deadline- 
first service in heavy-traffic acyclic networks. Ann. Appl. Probab. 14 1306-1352. 

MR2071425 

[24] Kruk, L., Lehoczky, J. P. and Shreve, S. (2006). Accuracy of state space coUapse 
for earliest-deadfine-first queues. Ann. Appl. Probab. 16 516-561. MR2244424 

[25] Kruk, L., Lehoczky, J. P., Ramanan, K. and Shreve, S. E. (2007). An explicit 
formula for the Skorokhod map on [0,a\. Ann. Probab. 35 1740-1768. 

[26] Kruk, L., Lehoczky, J. P., Ramanan, K. and Shreve, S. (2008). Double Sko- 
rokhod map and reneging real-time queues. In Markov Processes and Related 
Topics: A Festschrift for Thomas G. Kurtz. Inst. Math. Stat. Collect. 4 169- 
193. IMS, Beachwood, OH. MR2574231 

[27] Lehoczky, J. P. (1996). Real-time queueing theory. In Proc. of IEEE Real-Time 
Systems Symposium 186-195. IEEE Computer Society Press, Los Alamitos, CA. 

[28] LiMiC, V. (2000). On the behavior of LfFO preemptive resume queues in heavy traffic. 
Electron. Comm. Probab. 5 13-27 (electronic). MR1736721 

[29] Panwar, S. S. and Towsley, D. (1992). Optimality of the stochastic earliest dead- 
line policy for the G/M/c queue serving customers with deadlines. In Second 
ORSA Telecommunications Conference. ORSA (Operations Research Society of 
America), Baltimore, MD. 

[30] Prokhorov, Y. (1956). Convergence of random processes and limit theorems in 
probability theory. Theory Probab. Appl. 1 157-214. 



EDF QUEUES WITH RENEGING 



65 



[31] Ramanan, K. and Reiman, M. I. (2003). Fluid and heavy traffic diffusion lim- 
its for a generalized processor sharing model. Ann. Appl. Probab. 13 100-139. 
MR1951995 

[32] Ross, S. M. (1983). Stochastic Processes. Wiley, New York. MR683455 

[33] Ward, A. R. and Glynn, P. W. (2003). A diffusion approximation for a Markovian 
queue with reneging. Queueing Syst. 43 103-128. MR1957808 

[34] Ward, A. R. and Glynn, P. W. (2005). A diffusion approximation for a GI/GI/1 
queue with balking or reneging. Queueing Syst. 50 371-400. MR2172907 

[35] Whitt, W. (1971). Weak convergence theorems for priority queues: Preemptive- 
resume discipline. J. Appl. Probab. 8 74-94. MR0307389 

[36] Whitt, W. (2002). Stochastic- Process Limits: An Introduction to Stochastic-Process 
Limits and Their Application to Queues. Springer, New York. MR1876437 

[37] Zhang, J. and Zwart, B. (2008). Steady state approximations of limited processor 
sharing queues in heavy traffic. Queueing Syst. 60 227-246. MR2461617 



L. Kruk 

Department of Mathematics 

Maria Curie-Sklodowska University 

Lublin 

Poland 

AND 

Institute of Mathematics 
Polish Academy of Sciences 
Warsaw 
Poland 

E-MAIL ; Ikruk@hektor.umcs .lublin.pl 
K. Ramanan 

Division of Applied Mathematics 
Brown University 
Providence, Rhode Island 02912 
USA 

E-MAIL: Kavita_Ramanan@brown.edu 



J. Lehoczky 

Department of Statistics 
Carnegie Mellon University 
PiTTSBURCiH, Pennsylvania 15213 
USA 

E-MAIL: jpl@stat.cmu.edu 



S. Shreve 

Department of Mathematical Sciences 
Carnegie Mellon University 
Pittsburgh, Pennsylvania 15213 
USA 

E-MAIL: shreve@andrew.cmu.edu 



